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ABSTRACT 

The Orion-Eridanus superbubble is the prototypical superbubble due to its proximity and evolutionary state. 
Here, we provide a synthesis of recent observational data from WISE and Planck with archival data, allow¬ 
ing to draw a new and more complete picture on the history and evolution of the Orion-Eridanus region. We 
discuss the general morphological structures and observational characteristics of the superbubble, and derive 
quantitative properties of the gas- and dust inside Barnard’s Loop. We reveal that Barnard’s Loop is a complete 
bubble structure which, together with the A Ori region and other smaller-scale bubbles, expands within the 
Orion-Eridanus superbubble. We argue that the Orion-Eridanus superbubble is larger and more complex than 
previously thought, and that it can be viewed as a series of nested shells, superimposed along the line of sight. 
During the lifetime of the superbubble, HII region champagne flows and thermal evaporation of embedded 
clouds continuously mass-load the superbubble interior, while winds or supernovae from the Orion OB associ¬ 
ation rejuvenate the superbubble by sweeping up the material from the interior cavities in an episodic fashion, 
possibly triggering the formation of new stars that form shells of their own. The steady supply of material into 
the superbubble cavity implies that dust processing from interior supernova remnants is more efficient than pre¬ 
viously thought. The cycle of mass-loading, interior cleansing, and star formation repeats until the molecular 
reservoir is depleted or the clouds have been disrupted. While the nested shells come and go, the superbubble 
remains for tens of millions of years. 


1. INTRODUCTION 

The Orion-Eridanus superbubble is a nearby (~400 pc) ex¬ 
panding structure that is thought to span 20° x 45° on the sky 
( [Reynolds & Qgden|1979||Baily|2008||Pon et al.|2014a) . Be¬ 

cause of its proximity and evolutionary stage, it is traced at 
a multitude of wavelengths (Heiles et al. 1999) , serving as a 
benchmark to the study of superbubbles. The expansion of 
the bubble is most likely connected to the combined effects of 
ionizing UV radiation, stellar winds, and a sequence of super¬ 
nova explosions from the Orion OB association, Orion OBI 
<[bt aauw p%4[|Brown et al.|l994) . 

The line of sight expansion of the Orion-Eridanus super¬ 
bubble is previously estimated at ~ 15 km s 1 through line¬ 
splitting of Her, which also revealed a total ionized gas mass 
of 8 x 10 3 4 M q and a kinetic energy of E k in > L9 x 10 50 erg 
(Reynolds & Ogden 1979) . Large-scale spectroscopic map¬ 
ping of the entire Orion-Eridanus region in Hi reported a 
larger expansion velocity of the superbubble of ~ 40 km s' 1 
with a mass of 2.5 x_ 10 5 M 0 , containing a kinetic energy 


of 3.7 x 10 51 erg (Brown et al. 1995), which is shown to be 
consistent with the integrat ed mechanical lum inosity exerted 
by Orion OBI (~ 10 52 erg; Brown et~ak|l995 1 . The expan¬ 
sion of the superbubble is ailso traced through high-velocity 


superbubble (Burrows et; 

Here, we present observations of recent all-sky surveys with 
WISE and Planck, and combine them with existing sky sur¬ 
veys to provide a new and more complete insight in the his¬ 
tory and evolution of the Orion-Eridanus region. In particular, 
we will reveal that Barnard’s Loop is part of a complete bub- 
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ble structure, separate from the Orion-Eridanus superbubble, 
that sweeps up the mass-loaded interior of the pre-existing 
superbubble. It is argued that the Orion-Eridanus superbub¬ 
ble is larger and more complex than previously thought, and 
that the entire morphological appearance of the superbubble 
can be viewed as a series of nested shells, superimposed along 
the line of sight. The shells originate from explosive feedback 
from Orion OB 1 that accelerate, sweep-up, and compress the 
superbubble interior plasmas in an episodic fashion to form 
nested shells within the Orion-Eridanus superbubble. We ex¬ 
plore the origin of the shells, their relation with the subgroups 
of Orion OBI, and their impact on the molecular clouds and 
star formation efficiency within Orion. We discuss our find¬ 
ings in terms of the long-term evolution of the superbubble. 
We present the obervations in Sec. [2]and our results in Sec . [3] 
In Sec. [4j we discuss our findings, and summarize in Sec. [5^ 

2. OBSERVATIONS 

We made use of the all-sky surveys from Planck (|Planck| 
[Collaboration et al.|2014c), the Wide-Field Infrared Explorer 


(WISE; Wright et al. 2010]), the Leiden/Argentina/Bonn 


survey of Galactic H I "(LA"l3; |Kalberla et al.j [2005] ), 
Rontgensatellit soft X-ray background (ROSAT SXF 


and intermediate-velocity 

^as in several lines of sight towards 

Orion (Cowie et al. 

1979 

Welty et al. 2002). Soft X-rays 

emanate from the mi 

lion-degree plasma in the interior of the 


Snowden et al. |1997) , and an all-sky H er map that combines 

Finkbeiner 2003), including the 
;[Gaustad et al. 


several large-scale surveys 


Southern H-alpha Sky Survey Atlas 
2001), the Virgin ia-Tech Spectral-line Survey (VTSS; |Dei> 
nison et aI7||1998), and the Wisconsin Ha Mapper (WHAM; 
Haftfaer et al.|2003) . 


3. RESULTS 


Figure 1 reveals the large-scale structure of the Orion- 
Eridanus superbubble. In Her, the region exhibits various fila¬ 
mentary structures, including B arnard’s loop (Barnard|1894 i 

and the Eridanus filaments (e.g., Pon et al. 20 Projected 
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Fig. 1. — (a) Three-color image of the Orion-Eridanus superbubble. Ho- (blue) reveals ionized gas, surrounded by a layer of polycyclic aromatic hydrocarbons 
or stochastically-heated small grains, traced by the WISE 12 //in band (green). Red is Planck 353 GHz that probes columns of cold (20 K) dust. The polycyclic 
aromatic hydrocarbons/small grains trace the UV-illuminated edges of the clouds containing cold dust, resulting in that red and green blend to form yellow . To 
the north-east (top part the image) lies the Galactic plane, while to the north-west (right part of the image) lie high-latitude clouds (e.g., |Schlafly et al.|20i4j i. 
The numbered boxes define the apertures discussed in Sec. |3.2| while the light-green labelled solid lines mark the positions of cross-cuts used in Fig7(2] (b) 
Same region as panel (a), but only showing the Hu emission using a different stretc h to acce ntuate the faint shell structure which envelops the Orion region in 
Ho-. The white dot reveals the flux-weighted centre of the Ho bubble (Reynolds & Ogden 1 1979f (c) HI emission from the LAB survey, integrated in velocity 
space along the full extent of the Orion-Eridanus superbubble, -40 km s 1 < ty.sR < 40 km s 1 (Brown et al. 19951. Overlayed are contours of ROSAT 0.25 keV 
(orange dashed) and ROSAT 0.75 keV diffuse X-rays (light blue dashed). The 0.75 keV emission projected within Barnard’s Loop originates from bright sources 
within the ONC and IC 434 regions, but is blended with the more diffuse emission from the superbubble because of smoothing of the ROSAT maps. Below the 
white solid line the image has a different scaling to accentuate the faint HI filament intersecting apertures 2-4. Also overplotted is the shape of the superbubble 
aperture (S.A.; see text) (d) Schematic representation of the region. Her gas structures are shown in blue lines/contours, while dust structures are plotted in red. 
The solid lines trace faint filaments of the region to guide the eye in panels (a) & (b). The stars that form the well-known constellation of Orion a re also p lotted to 
appreciate the immense size of the region (black diamonds). Members of the different subgroups of the Orion OB association (Blaauw|l964||Brown et al.| i 994) 
with spectral type B2 or earlier are also shown in different grey symbols shown by the legend. See text for the explanation of the different abbreviations. 
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on top of Barnard’s Loop are the Orion A and Orion B molec¬ 
ular clouds (OMC A and OMC B), the Orion Nebular cluster 
(ONC) associated with the Orion nebula, and the IC 434 emis¬ 
sion nebula that is characterized b y a cha mpagne flow of ion¬ 
ized gas (|Tenorio-Tagle| 1979 1 Sec. |3.3.4)». Another prominent 
feature is the A Orionis bubble, a 10° circular Hu region, en¬ 


capsulated by a swept-up shell of gas and dust. Below, we will 
refer to different parts of the superbubble in RA/Dec space 
(see the orientation in Fig. [T]). 

Figure [2] shows cross-cuts through the Orion-Eridanus su¬ 
perbubble, revealing the gas- and dust structures of the super¬ 
bubble and its interior components. Note that Barnard’s Loop 
peaks dramatically in the ionized gas, while in the dust and 
neutral hydrogen it is only seen as a modest increase in the 
Planck, WISE, and LAB channels (Fig. [2]i,b,c). At this point, 
we find no evidence that a dense shell encapsulates the ionized 
portion of Barnard’s Loop, which presence was predicted by 
Cowie et alT|( |1979| (Sec. [4}. In contrast, the cross-cut shown 
in Fig. [2jl that traces the western bubble shell at high latitudes 
reveals a stratified structure, where a dust/ HI shell is clearly 
observed to encapsulate the Eridanus C filament. The absence 
of an HI shell associated with Barnard’s Loop is surprising, 
and encouraged us to investigate the gas- and dust content of 
Barnard’s Loop and the Orion-Eridanus superbubble in more 
detail; our findings are presented in the next subsections. 


3 . 1 . Dust in Barnard’s Loop 
3 . 1 . 1 . Optical depth 


Heiles et al. (2000) defined two small regions of Barnard’s 
Loop which are clear of any dense molecular clouds along 
the line of sight; we use the same apertures in the following 
analysis (’BL top ’ and ’BLbot’; Fig. HTl to characterize the dust 
contained in Barnard’s Loop. The distance toward Barnard’s 
Loop is uncertain; here, we will assume that B arn aid/s_ Loop 
lies at the distance of the Orion nebula (400 pc; |Sandstrom 
et al.||2007] |Menten et al.||2007] >. Electron de nsities inside 
Barnard’s Loop are estimated at n e - 2.0 cm -3 (Heiles et al. 


2000 1 , and at n e — 3.2 cm " 3 f rom the surface brightness of 
the H/f fine ( |Q’Dell et al.|2011~| . A similar analysis of the Her 
surface brightness here leads to n e - 3.4 cm~ 3 , after cor rect- 
ing for limb-brightening (a factor 7; O’Dell et al.|2011 1 of a 
constant density shell with an apparent radial thickness of 1 °, 
which amounts to 7 pc at the distance of 400 pc. 

Using n e = 3.4 cm ~ 3 and a radial thickness of 7 pc, adopt¬ 
ing a normal gas-to-dust ratio of Ah/A v = 1.9 x 10 21 cm ~ 2 


magnitude (Bohlin et al. 1978), where Ah is the column 


density of hydrogen and A v the visual extinction, and a UV- 


to-visual dust absorption ratio of 1.8 (Tielens & Hollenbach 
1985), the radial dust optical depth for UV photons with en¬ 


ergies between 6 and 13.6 eV is Tjjv = 0.07. 

Through a modified-blackbody fit to the Planck data, the 
Planck R1.20 dust model provides all-sky maps of dust opti¬ 
cal depth (at 353 GHz), temperature, spectral index, and the 
total dust emission integrated over frequency. We obtain an 
observed optical depth of Barnard’s Loop at 353 GHz, T 353 > 0 b s , 
by correcting the value of the optical depth from the Planck 
model at the location of Barnard’s Loop for a background 
value. For this method we use the bins defined in Figs. [2]i&c: 
as Barnar d’s L oop is not encapsulated by a dense sheIl(Fig. 
[2]and Sec. |3.2[ >, we determine the background in a region just 
outside of the ionized portion of Barnard’s Loop (which is 
significant; ~ 60% of the total value for both regions). In this 
way, we isolate the dust grains contained within the ionized 


part of Barnard’s Loop. Then, the observed optical depths 
are T 35 3i0bs = 9.0 x 10 ~ 6 (for BL top ) and r 3 53 >obs = 5.1 x 10 “ 6 
(BLbot)- We can compare this to the expected optical depth, 
053,cal, from the density (3.4 cm -3 ), estimated line of sight 
depth (45 pc for a 1° thickness) and the opacity at 353 GHz 
(cr v j= 7.9 x IQ - 27 cm 2 H 1 (appropriate for the diffuse ISM 
(Planck Collaboration et al. 2014a| l. We calculate 7353 ^ = 
4.0 x 10~ 6 , in close agreement with T 353 j 0 b s , which implies 
that dust mixed with the ionized shell of Barnard’s Loop is 
alike to dust seen in the diffuse ISM. 


3 . 1 . 2 . Temperature 


Heiles et al.] ( |2000| ) concluded that grain temperatures in 
Barnard’s Loop are higher than in HI regions, and proposed 
that Lya photons might contribute to the heating of the dust. 
However, direct heating from the radiation field of Orion OB 1 
was not considered in their study. We derive an intensity of 
Go = 14 reaching Barnard’s Loop in units of the Habing field 
( |Habingj)l968| , taking into account the radiation from the 7 
most massive stars in Orion OB (9 Ori C, 6 Ori, t Ori, 6 Ori 


A, cr Ori, e Ori, C Ori), using stellar parameters from Martins 
|et al.| ( |2()()5[ >. However, |Q’Dell| ( f2001) l shows that radiation 
from the Orion nebula (6 Ori C and 6 Ori A) might not reach 
Barnard’s Loop, as the nebula is optically thick to ionizing 
radiation in all directions, except possibly to the southwest. 
This would decrease Go to 11. 

The observed dust temperature from the Planck model in 
BL top and BLbot are Tjj = 18.7 K and 7d = 20.5 K. These val¬ 


ues are consistent with temperatures found by Heiles et ah] 
(2000 1 , 7d = 19.6 K and 7/j = 20.3 K for BL top and BLb ot , re¬ 
spectively, and the mean value 7/j = 19.7 K for the entire sky 
(Planck Collaboration et al.j |2014a l, even though Go in BL 
insignificantly raised above the average interstellar radiation 


field of the solar neighborhood (= 1.7 Habing field; Tielens 
2005]). For comparison, at Go = 14 a ‘typical’ 0.1 pm silicate 
(graphite) grain would be heated towards 7d = 19.2 (22.5) K 
(Tielen s]2005] >. At Gq_= 11, this would be 7 d =18.6 (2L8) K. 


In this respect, Planck Collaboration et al. (201 4a)> and Planck 


Collaboration et al.|(|2014b|" argue that observed dust temper¬ 
atures do not simply trace the intensity of the radiation field, 
but in addition reflects changes in grain properties, such as the 
size distribution, grain structure, and material changes. In any 
case, for both BL top and BLbot, we do not see the necessity 
of an extra heating source to account for the observed grain 
temperatures, such as Lya heating (Heiles et al. [2000| . 


3 . 2 . Dust and gas in the Orion-Eridanus superbubble 
3 . 2 . 1 . Global observational characteristics 

The all-sky surveys listed in Sec. fallow us to determine 
the global observational characteristics of the Orion-Eridanus 
superbubble, which is of potential use to the study of more 
distant superbubbles that lack the spatial information of the 
Orion-Eridanus region. We use apertures that encompass the 
entire superbubble, the OMCs, and Barnard’s Loop, respec¬ 
tively. For the OMCs and Barnard’s Loop, these apertures 
coincide with the outlines depicted in Fig. [III. While the 
full extent of the Orion-Eridanus superbubble will be thor¬ 
oughly discussed in the remainder of this paper, here we al¬ 
ready choose the ‘entire superbubble aperture’ to encompass 
the majority of the Ha emission observed in Fig. |TJ>, while 
tracing the outside of HI filaments seen at -40 km s' > vlsr 
> 40 km s 1 (Fig. fTp), corresponding to the velocity extent 
in which the Orion-Eridanus superbubble has previously been 
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Fig. 2.— Profiles of cross-cuts defined in Fig |TJi. Cuts start at the labelled side in Fig. [TJi. (al) Along the cut from left to right: Barnard’s loop; the Orion 
B molecular cloud (OMC B) and the IC 434 emission nebula, the GS206-17+13 shell; tne excavated bubble interior; and the MBM 20 high-latitude cloud, 
overlapping the Eri A Hn filament. The dashed vertical lines highlight several prominent features seen in the gas (blue) and in dust (red). (a2) Inset of the IC 434 
emission nebula Hn profile, plotted as intensity (in units o f Rayl eigh) as a function of distance from the ionizing star r (in degrees). The emission profile exhibits 
the structure of a champagne flow of ioniz ed gas {Tenorio-Tagle11979). We have marked the location of the ionizing star, <x Ori (r*), the shock front ( r s ), and 
the rarefaction wave (r rw ). See Sec. |3.3.4| for a full description or the champagne flow structure, (b & c) Cuts through Barnards’s loop (BL top and BLbot; |Heiles| 
|et al.|2000), revealing a single peak most p rominently seen in H a. The dashe d lines on the bottom define the bins where the mean optical depth at 353 GHz, 
r 353 ,obsi from the Planck dust model R1.20 {Planck Collaboration et al.|2014a) is measured (see text): ‘BL’ is the bin for Barnard’s Loop, ‘BG’ is the bin for the 
background determination, (d) Cut through tne western edge ot the superbubble, coincident with the Eri C filament. Here, the gas and dust reveal a stratified 
structure, where the ionized gas peaks inside the dust shell. 


identified ( |Brown et al.|1995] >. We measure the Her luminos¬ 
ity and total far-IR to sub-mm luminosity contained within 
the apertures, which are compared with stellar parameters of 
Orion OB 1. 

For the entire superbubble, we proceed as follows. The 
north- and east side of the superbubble aperture extend close 
to the Galactic plane and high-latitude clouds. The encom¬ 
passed IR emission will therefore be contaminated by emis¬ 
sion along the line of sight not related to the superbubble. 
While disentangling the separate contributions to the IR is not 
straightforward, we attempt to isolate the emission from our 
region of interest by measuring the average far-IR integrated 
intensity, /irhi, within the aperture labeled as ‘IR fil.’ (Fig. 
[TJl), i.e., the dusty shell that encapsulates the limb-brightened 
Eridanus B and C filaments known to be associated with the 
superbubble. The intensities are taken from the Planck R1.20 
dust model ( |Planck Collaboration et al-120 14a| ) that uses IRAS 
100 /rm and Planck 857, 545, and 353 GHz data to determine 
the total integrated far-IR intensity. Subsequently, we define 
emission from regions that are much brighter than observed 
in the limb-brightened filaments, i.e., 7 ir > 3/|r.i,i, as not be¬ 
ing part of the superbubble. This procedure provides a mask 
tracing the regions of high IR brightness towards the Galac¬ 
tic plane, and high-latitude fore- and background clouds well. 
Note that this procedure filters out the emission from, e.g., 
the OMCs as well; we ultimately add the contributions of the 
OMC and Barnard’s Loop apertures to get the full IR inten¬ 
sity of the superbubble and the relative contributions of the 
components (see Tab. [TV For the Ha map, the background 
contamination is muchless severe and a masking procedure 
is not necessary. In this respect, we will argue below that all 
of the observed Ha in Fig. [7]i is caused by photo-ionization 
from Orion OB 1. 

While the masking procedure accounts for the high-IR in¬ 
tensities of contaminating clouds towards the Galactic plane, 
it does not correct for the diffuse low-level background ob¬ 
served over the entire images. Therefore, we apply a global 
background subtraction to both the far-IR and Ha maps by 
measuring the background level within a circular region of 
0.4 degree radius at (l,b) = (229.6, -33.7) outside of the super¬ 
bubble aperture that appears to be free of obvious emission 
features and provides a proper representation of the diffuse 
background emission averaged over latitude. For both maps. 


we do not include the A Ori region in our analysis as it appears 
to be encircled by an ionization- bounded H II region (Warren 
& HesserJT977 Wall et al.j 1996 >, and can be regarded as be- 
ing separated from Orion OB 1. 


ing separated] 

Table [I] lists two different values of Lir. The first value 
is obtained by filling the masked region with the average IR 
brightness from the unmasked part of the superbubble aper¬ 
ture. In this way, we get an estimate of the full IR intensity of 
the entire superbubble if it were unaffected by fore- or back¬ 
ground emission along the line of sight. The second value 
is obtained by setting the masked part to zero, and therefore 
neglects the masked part of the superbubble. The choice of 
including or excluding the masked part of the superbubble af¬ 
fects the inferred luminosities from the OMC and Barnard’s 
Loop as well, as these apertures are completely immersed in 
the masked region of the image, and the chosen intensity of 
the masked region thus effectively acts as a ‘global’ back¬ 
ground for both regions. Below, we will include the masked 
region in our interpretations, as we believe that the inclusion 
of this part of the superbubble is essential to the analysis of the 
total energetics of the region. However, we refer the reader to 
Table [TJ for the implications of this choice and its impact on 
the derived parameters. The total far-IR luminosity is Lir = 
4-7:I tot S , where 7 tot and S are the total integrated intensity and 
projected surface area of the emitting region. For simplicity, 
we assume that the entire bubble is located at 400 pc. The 
measured IR luminosity is compared with stellar parameters 
of Orion OBI through L* = L IR /L*, where L* is the total lu¬ 
minosity of the Orion OB region (1.7 x 1 0 6 L 0 , re-evaluated 


at d — 400 pc, an d excluding A Ori; Warren & Hesser 1977 
Wall et al. 1996[>. The parameter defines the fraction of 


stellar radiation that is captured by the dust and re-emitted in 
the IR averaged over all solid angles. 

The Lya photon rate, N\_ y „, is calculated by converting the 
observed Ha photon rate to a Lya photon rate using the ratio 
o'b/o'Hq-' where a Ha is the effective recombination coefficient 
of Ha (1.31 x 10~ 13 cm 3 s _1 at T - 6000 K, which is the 
electron temper ature in Barnard ’ s Loop and the Eridanus fil- 
aments from |Heiles et al.||2000] IQ’Dell et al.||20H~l |Madsen| 
|et al.|[2006l we use this value throughout the superbubble), 
and a B islhe total recombination coefficient of hydrogen to all 
levels but the ground state (2.6 x 10 


Osterbrock 


& Ferland 2006 1 . Case B recombination requires that N\ 


Ly a 
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TABLE 1 

Global observational characteristics of the Orion-Eridanus superbubble 


Region 

T A 

(K) 

Lir 

(10 5 Lo) 

?* 

L]_ya 
(10 5 Lo) 

£ion 

Superbubble 

19.5 

7.9 (5.7) (a) 

0.47 (0.34) 

1.08 

0.94 

OMC 

17.9 

2.8 (2.9) 

0.16(0.17) 

0.16 

0.14 

Barnard’s Loop 

19.6 

0.5 (0.6) 

0.03 (0.04) 

0.20 

0.17 


Note. — Listed are: the dust temperature, T&\ the infrared luminosity, Lir; 
fraction of Lir to the total luminosity of Orion OB 1, ; luminosity measured 

from Ho and converted to Lyo, L\ ya ; fraction of Lya photons to the total 
amount of ionizing photons of Orion OBI, trion- The values between brackets 
denote the values measured when the masked part of the superbubble is set 
to zero (see text). ( ' L) : The total IR luminosity in the superbubble + OMC + 
Barnard’s Loop apertures. 


equals the total amount of recombinations, which is a quan¬ 
tity that can directly be compared to Orion OB 1 through £)<,„ 
= Mya/Wion, where Mon is the total number of ionizing pho¬ 
tons from Orion OBI (2.7 x 10 49 ph s _1 , O’Dell et ah 2011). 
This ratio measures the fraction of ionizing photons captured 
by the gas and converted to Lya. In order to directly compare 
with the IR luminosity Lir, we convert Nh ya to a luminosity 
through Liya = Ni ya hv a , where hv a is the Lya photon en¬ 
ergy, and define L ion = N ion hv a , such that £ on = N Ly JN ion = 
L[ yfylL um . 

The calculated luminosities, Lir and Li ya , are denoted in 
Tab. [I] Table Q] shows that half of the total amount of stel¬ 
lar radiation of Orion OBI is trapped in the superbubble and 
re-radiated in the IR, ~ 0.5. The molecular clouds account 
for about 35% of the total IR emission of the superbubble; the 
contribution from Barnard’s Loop is negligible. In contrast, 
£i on of the OMC (0.14) and Barnard’s Loop (0.17) are roughly 
similar, which is because the Orion nebula and IC 434 emis¬ 
sion nebu la are currently breaking out of the molecular clouds 
(Sec. |3.3| > and are contained within the OMC aperture. Even¬ 
tually, all ionizing photons are absorbed within the superbub¬ 
ble aperture (L mn ~ 1). This implies that, on average, Orion 
OB 1 can provide the necessary ionizing power to illuminate 
the Ho- structures detected in Fig [lj:, including the Ha fila¬ 
ments that run along the outer edge of the superbubble aper¬ 
ture. Nonetheless, it is possible that individual small-scale 
stmctures may still be too bright given their size and distance 
from Orion OBI ( |Pon et al.|2014b| l. Note that some of these 
filaments lie outside of what has previously been thought to 


be the edges of the Orion-Eridanus superbubble Reynolds 
& Ogden|[l979l |Heiles et al.|fl999| [Bally|[20081 lfe_ 


on et al. 

2014a). The dust temperatures measured inside the superbub 


ble aperture and Barnard’s Loop are similar to that observed 
for th e entire diffuse sky by Planck (see discussion in Sec. 
|3.1.2|i. Inside the OMC, this value is somewhat lower, which 


is what is expected for dense regions (Planck Collaboration 
|et al.|2014al >. 


3 . 2 . 2 . Tracing the superbubble structure through dust and gas 

To investigate the gas- and dust content in specific regions 
throughout the Orion-Eridanus region, we use apertures de¬ 
picted in Fig. [l]that are defined such that the majority of Ha 
and far-IR emission from the regions are enclosed (for exam¬ 
ple, the stratified emission from the limb-brightened emission 
of the superbubble wall; Fig. [2]f), while a low background 
level is ensured through comparison with the Planck and LAB 
maps (the Ha emission does not show a high background in 
general). As the diffuse background emission varies signifi¬ 
cantly over the entire Orion-Eridanus region, here we use lo¬ 


cal background values to properly measure the luminosities 
contained in the apertures scattered over the Orion-Eridanus 
region. For the apertures at b > -30 degrees, we define a back¬ 
ground at ( l,b ) = (233.2, -28.6) within a circle of 0.4 degree 
radius. For the apertures at b < -30 degrees, we define a back¬ 
ground at (Lb) = (221.9, -49.5). Region 5 is an exception to 
this rule, as the previously defined background levels were 
not representative of the local value because of its location 
amongst high-latitute clouds. For region 5, ( l,b ) = (176.3, - 
15.0) was chosen. 

Lya photons will be resonantly scattered many times in the 
HII region because of its large line cross-section, but even¬ 
tually be absorbed by dust, contributing to the heating of the 
dust. The ratio of heating rate by stellar photons, Tyv, to Lya 
photons, TLyo., can be written as (Tielen s]2005 i: 


I\jv _ 7ra 2 /7 d (2 a b s L* 

r Lya , 4nr 2 n 2 aBhv a 

Here, iter and n A are the geometrical cross-section and num¬ 
ber density of the grains, gabs is the average radiation ab¬ 
sorption efficiency of the dust, L* is the total luminosity of 
Orion OBI, r the distance to the source, and n is the hydrogen 
number density of the gas. Assuming equilibrium between 
photo-ionization and recombination, for a shell surrounding 
an empty cavity we have f on N lon = Anr 2 n 2 a^,Ar , where / on 
is the fraction of incident ionizing photons available that are 
absorbed by the gas locally. Then, Eq. [Tjreduces to 


T UV _ 7ta 2 n d Q abs Ar | 

( L * ] 

i - Td l 

( u \ 

r Lyor fon 


' /.on 1 

\N[ on hv a J 


( 2 ) 


where r d is the (radial) absorption optical depth of the shell. 
With L* = 1.7 x 10 6 L 0 and Mon = 2.7 x 10 49 ph s _1 (Sec. 
|3.2. 1| >, we have (L*/A lon /iv„) *15 for Orion OBI, and there¬ 
fore the ratio Tuv/TLyo- equals unity for T d /f on - 0.07. 

We can make a priori estimates of f on for the regions de¬ 
fined in Fig. [T] For Barnard’s Loop, we have measured an 
average £j on = 0.17 over the entire structure (Table [lj. A 
factor of 0.50 would be expected for a half-sphere geometry 
if Barnard’s Loop were optically thick to ionizing photons. 
Possibly, a fraction of the ionizing photons from Orion OB 1 
is trapped within the OMCs before reaching Barnard’s Loop, 
which would lead to Li„n.BL = 0.50 - Li 0n ,OMC = 0.36 (Table 
[TJ. Still, these numbers reveal that 50% - 66% of the ionizing 
photons pass through the half sphere encompassing Barnard’s 
Loop. Hence, we estimate that / on is in the range -0.33 - 
0.50 for the Barnard’s Loop apertures BLbot and BL top . The 
fact that Barnard’s Loop is optically thin for the ionizing flux 
of Orion OBI can also be inferred from the morphological ap¬ 
pearance of the cometary clouds L1617 and L1622, which are 
located behind Barnard’s Loop as measured from the Orion 
OB association and have sharp I-fronts on their sides facing 
the Belt stars (Bally et al. |2009| l. We note that the internal 
structure of Barnard’s Loop may be more complex than that 
of a homogeneously distributed half-sphere assumed here. A 
patchy structure containing optically thick clumps ( f OD - 1) 
and holes through which the photons would leak unhindered 
(/ion ~ 0) could, in principle, lead on average to an f on of or¬ 
der 0.33 - 0.50. However, Barnard’s Loop, similar to region 
1 and region 6, is not associated with HI. The remaining re¬ 
gions are associated with neutral hydrogen, where / on will be 
of order ~ 1 because of the small mean free path of ionizing 
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photons through Hi (Tielens|2005|. 

We measure both Lir and L\ yll within the apertures shown 
in Fig. [T] and compare these measurements with luminosities 
from Orion OBI. The structures contained in each defined 
aperture will subtend a solid angle Q as measured from Orion 
OBI. Now, any ionizing photon that is absorbed by hydrogen 
and converted to Ly a will eventually contribute to the heating 
of the dust and, therefore, the total IR luminosity that em¬ 
anates from each region can be written as 

LlR = Ojv + 0.y<T — T ( | L* G[iy + ./ionNj on /zV (f .O[ Vff ,, (3) 

where L uv is the dust luminosity provided by direct absorp¬ 
tion of stellar photons. In the classical picture of an expanding 
superbubble, a dense, neutral shell will wrap around the 10 4 
K region created by the reverse shock traced by bright Ha 
emission ( [Weaver et al.|1977] >. Therefore, we can assume that 
the IR and Ha trace structures which measure the same solid 
angle with respect to Orion OBI, Ojjv = C Lya ,, and we can 
write 


Oir/LQ 
(f^Ly al L ion) 


i"d / Nj ()n /iy n 

J7n + \ Li, 


(4) 


which describes the ratio between the amount of energy that 
is captured by the dust and re-emitted in the IR, to the amount 
of energy that is captured by the gas and converted to Lya. 

Our results are shown in Fig. [3] All studied regions show 
Lir > Li, ya (Fig. [Tji), which reveals that a (substantial) amount 
of stellar radiation is directly absorbed by the dust. Here we 
note that Lir = 2L l y(f when r ( j V = T| y „. (cf. Eq. |3|. An inter¬ 
esting trend can be seen from Fig. [3ji. The regions 2-5, po¬ 
sitioned on filamentary structures seen in gas and dust along 
the outside of the superbubble aperture depicted in Fig. T] re¬ 
veal significant higher fractions of Tjf on compared to the re¬ 
gions located in Barnard’s Loop (BLt, ot , BL top ). T his im plies 
that stellar photons, provided by Orion OBI (Sec. |3.2.1|), are 
more efficiently absorbed in regions 2-5 compared to that 
in Barnard’s Loop. In other words, the dust optical depth for 
stellar photons associated with Barnard’s Loop is very small. 
In Fig.]3j}, we rewrite Eq. [4]for T d // mn and, together with our 
estimations for f on as described above, plot the derived dust 
optical depth in the regions as a function of projected distance 
ufproj from (l,b) - (206.5,-18.0), the flux-weighted center of 
the bubble as measured f rom the Ha emission (see Lig. [IB; 
[Reynolds & Ogden|1979| ). We exclude A Ori as it seems tobe 
encircled by an ionization-bounded Hu region. Still, in Sec. 


3.3.3 we will show that the ionized gas in the A Ori region 
shows the characteristics of a champagne flow (Tenorio-Tagle 


1979)>, in which case the Hn region will be density bounded 
on the side of the outward champagne flow, and radiation will 
likely be able to escape the region. However, the circular mor¬ 
phology of the A Ori bubble detected in Ha (Fig. |4]» reveals 
that the cham pagne phase, if present, may not have fully de¬ 
veloped (Sec. |3.3.3[) and, therefore, we will assume that the A 
Ori HII region is still ionization-bounded. 

For the BL apertures, we plot t he results for T d versus c/ prl) | 
assuming f on = 0.5 (Sec. 3.2.2[ ). For regions 1 and 6, the 
absence of a clear neutral HI filament prohibited us from pre¬ 
dicting /; on , and only an upper limit to r,| is given. The results 
plotted in Fig. [3b, albeit suffering from low-number statistics, 
indicate a trendwhere Td rises with increasing </ )m j. 

From Fig. [3] it becomes apparent that dust associated with 
Barnard’s loop is not able to capture all of the radiation emit¬ 


ted by the Orion OBI association. In this context, the bright 
Ha filament that surrounds Barnard’s Loop, from region 1 
up to the Galactic plane (Fig. jT]», is lit up by ionizing pho¬ 
tons breaking through Barnard’s Loop. A straight line from 
the flux-weighted centre of the superbubble (Fig. [T}, through 
Barnard’s Loop to the outer filament in region 1, reveals a flux 
ratio of /-/i /L rC g.i ~ 5, where L/l and L reg | are the measured 
peak Ha fluxes in Barnard’s Loop and the outer filament in 
region 1, respectively. At a distance ratio of t/urM-cg.i ~ 2.5, 
this is consistent with both the structures absorbing a simi¬ 
lar amount of ionizing photons, confirming / on ~ 0.5 for the 
Barnard’s Loop apertures as was assumed in our derivations 
of Td above. 

As a sanity check to our methodology, we also plot in Lig. 
[3] the results for region 6 (corresponding to the southern end 
of the Eri A filament), which reveals low Td compared to the 
other Eridanus filaments. However, the association of the 
Eridanus A filament with the superbubble is still controver¬ 
sial and is discussed thoroughly in Pon et alT| (|2014b)>. For 
completeness, we show the results for apertures located on 
the G203-37 and GS204-31 clouds (region 7 & 8) which are 
thought to be inside and at the near-side of the Orion-Eridanus 
superbubble (Snowden et al.||1995|>. The clouds are highly 
opaque to 0.25 keV radiation wi th intervening hydrog en col¬ 
umn densities of ~8 x 10 20 cm 2 (Snowden et al. 1995), trans¬ 


lating to an optical depth a t UV w avelength of Tuv ~ 0.8 using 
the relations given in Sec. 3.1.1 in rough agreement with our 
measurements shown in Fig. |5] 

With /i on = 0.5 for the Barnard’s Loop apertures, the dust 
optical depth as measured from the IR (Td = 0.06 and Td = 
0.08 for BL top and BLbot, respectively; Fig. S 1 & b) is in 
agreement with a homogeneous shell containin g du st similar 
to that seen in the diffuse ISM (tuv = 0.07; Sec. |3.1| >. Finally, 
we note that while Lya contributes about 25-50% of the total 
energy absorbed by dust in Barnard’s Loop (Fig. |3|i), th is will 
have little effect on the dust temperature (Sec. [3. 1 ,2|» as 7 d 
only weakly depends on the heating rate. In particular, 7/ °c 
(Tuv + I/ycr) 1 ^ 6 for silicate grains (Tielens,2005). 


3.3. Expanding shells within the Orion-Eridanus region 
3.3.1. High- and intermediate velocity gas 

The presence of high-velocity (HV) absorption features of 
ionized gas between -120 km s _1 < vlsr -80 km fT 1 was 


detected in mu ltiple lines of sight towards Orion ( Cowie et al.| 
|1979| |Welty et al.||2002|). In addition, intermediate-velocity 


20 km s 1 ( 

Cowie & York 1978 Cowie et al. 1979 Huang 

et al. 1995; 

Welty et al .75002). Both components show sim- 

ilar depletions, ion ratios, and physical conditions, albeit the 
IV gas contains column densities that exceed that of the HV 

gas by a fac 
dubbed by 

itor of 3 (e.g.. 

Welty et al. 2002). The HV gas. 

Cowie et al. C 

979) as “Orion’s Cloak”, has an 


angular diameter of at least 15° and is ascribed to a recent su¬ 
pernova that, from its dynamical expansion scale, should have 
occurred some 3 x 10 5 yr ago. 

Figure |4] shows Ha velocity map s from the WHAM Ha 
spectral survey (Haffner et al. 2003[ > of the Orion-Eridanus 


superbubble, covering the velocity space between -80 km s 
< i' LSR < +80 km s 1 . Overplotted are sightlines fro m ab - 


sorption line studies ( 

Cowie & York 1978 

Cowie et al. 

1979 

Huang et al.]jT995 V 

Velty et al. 2002). The HV and! 

V gas 


extend outside of the traditional limits of the Orion-Eridanus 
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Fig. 3.— Apertures tracing the structure of the Orion-Eridanus superbubble, (a) Infared versus Ly a luminosities measured in the apertures shown in Fig. IT] 
The solid line shows where the infrared and Ly a luminosities are equal. The lower dashed line shows where the infrared and Ly a heating rates are equal, 
corresponding to rJf on = 0.07 (see text), where us the dust optical depth at UV wavelengths and f[ on is the fraction of ionizing photons available absorbed 
by the gas. The dotted lines are similar, but for increasing rd. Overplotted are the luminosities measured in the apertures drawn in Fig. [T] color-coded as follows: 
regions 1-5 ( red) are positioned along the projected outer shell of the superbubble aperture (Fig. 171:); Barnard’s loop (BL top and BLbot) ( blue ); region 6 - and 
8 (gray) are the Eridanus A filament and molecular clouds close to the Orion OB association, (bji he radial dust optical depth at UV wavelengths versus the 
projected distance towards the Ht* flux-weighted center of the superbubble. 


superbubble ([Reynolds & Ogden||1979| 

Bally|2008| Pon et al |2014a| ), as demon 


Heiles et al. 1999 
emonstrated by the sight- 
lines of k Ori, A Ori, and A Fri. While it must be recognized 
that the distribution of the gases are somewhat patchy (Cowie 
et al. 1979| >, it is clear that the IV and HV gas trace two dls^ 
tinct components. Moreover, the HV gas seems to be confined 
to a more limited area compared to the IV gas, as reflected by 
the sightlines that show only the IV gas. 

We note upfront that electron temperatures within 
Barnard’s Loop and A Ori are determined at ~6000 - 7500 


K 


eiles et al. |2000l IQ’Dell et al.|201 H|MacL 


sen et al.| |2006| > which excludes that the HV gas originates 
from high-velocity wings of a thermal Gaussian line profile 
(the thermal broadening can be estimated with y/2kT/mu ~ 
10 km s" 1 , where k is Boltzmann’s constant and ;;jh is the 
mass of a hydrogen atom). The velocity information of the 
Eridanus filaments are discussed in |Pon et al7| ( |20 14b| >. Of 
particular interest here are the encircled regions towards the 
east of the superbubble, covering the Barnard’s Loop region 
and the A Ori region, as well as smaller-scale bubbles. 

3 . 3 . 2 . Barnard's Loop bubble 

The intensity peak of Barnard’s Loop lies between -20 km 
s 1 ^ 1 'lsr ^ +20 km s -1 , but the structure is visible through¬ 
out the entire velocity space covered by the WHAM survey. 
Here, we reveal the presence of an ionized filament, visi¬ 
ble only at positive velocities vlsr > +20 km s 1 . This fil¬ 
ament appears to form a complete bubble structure together 
with the bright crescent of Barnard’s Loop, centered on ( l,b) 
= (193,-20) with radius 14°. The lines of sight covered by 
the absorption survey of Cowie et al.|(|l979) indicate that the 
HV gas originates from the encircled Her bubble structure, 
or the ‘Barnard’s Loop Bubble’ (Fig. |4|), suggesting a con- 


mated the size of Orion’s Cloak at a minimum of 15 degrees 
that matches the size of the Barnard’s Loop bubble. Thus, 
the Barnard’s Loop bubble may well be part of the supernova 
remnant associated with Orion’s Cloak. 

The most pronounced radial velocity effects are expected 
to originate from the center of an expanding bubble. For the 
HV gas, these velocities lie outside of the velocity range cov¬ 
ered by the WHAM survey (-80 km s 1 < vlsr £ 80 km 


Haffner et al. 


2003), rendering WHAM unable to pick 


up the H V gas detected by C opernicus (Cowie et al.||1979]> 


and HST (Welty et al. 2002), although it is questionable if 
WHAM would be able to pick up the low surface bri ghtness 
of the HV gas expected in Ha (~3 x 10 3 Rayleigh; 


Cowie 


et al.|| 1979] >. The OMCs are roughly centered on vlsr = 10 


km s 1 in velocity-space (Wilson et al.||2005i; if we assume 


a systemic velocity towards the Orion region of vlsr = 10 
km s 1 and connect the HV gas with the expansion of the 
Barnard’s Loop bubble, we estimate its expansion velocity at 
1 . However, it is not clear how much of the 


t'exp = 100 km s 


mass within the Barnard’s Loop bubble (Sec. 4.3.1 1 is associ¬ 
ated with this expansion velocity. The coverage of sightlines 
of the Orion-Eridanus region is limited, and the HV gas could 
in principle trace small, high-velocity components related to 
the Barnard’s Loop bubble, rather than tracing the bulk ex¬ 
pansion of the bubble itself. Here we implicitly assume that 
the bulk of the gas in the expanding shell is expanding at the 
adopted expansion velocity. The pre-shock densities for the 


Welty et al. 2002) 


HV gas towards C, Ori (no = 3 x 10 J cm 
then implies that the bubble expands either into the 
tercloud Medium (HIM; T ~ 10 6 K, n ~ 3 X 10 3 cm 


I 


nection between both components. Cowie et al. (1979) esti 


lot In 
~ 3 ) that 

pervades the Galactic Halo, or a pre-existing cavity, such as 
the superbubble interior. We defer a thorough discussion of 
the Barnard’s Loop bubble and its connection to the HV gas 
to section [4] 
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Fig. 4. — WHAM Her observations of the Orion-Eridanus superbubble. Shown are color-coded intensity maps (in units of Rayleigh) integrated over a velocity 
range in vi ,sr that is indicated on top of each panel. Note the different scaling of each individual panel. Overplotted are lines of sight from absorption line studies 
(Cowie & York 1978, Cowie et al. 1979, Huang et al. 1995, Welty et al. 2002 1 that show both high-velocity and intermediate-velocity gas (inverted triangles), 
lines of sight that show only intermediate-velocity gas (diamonds), and lines of sight that do not show a clear detection of high-velocity or intermediate-velocity 
gas (squares). The classes are plotted in different gray scalings to make them discernable, and we have explicitely labeled some of the stars in the lower right 
panel for orientation. The stars 7r 5 Ori and e Ori are denoted in italics as observations onl y covered ve locities of |vlsrI < 60 km s 1 and |i>lsrI < 100 km s 1 . 
respe ctively ( jCowic et al. 119791. Measurements of HV/IV gas in A Lep are hindered {Cowie et al.| 1979) because of stellar rotation (s), while y Ori lies at 77 pc 
I van I.ceuwen|2007) and may trace foreground gas (/). 


3 . 3 . 3 . A Ori bubble 

The circular-symmetric H II region surrounding A Ori has 
been known for a long time. Throug h a st udy of the star for¬ 
mation history of the A Ori region, |Dolan & Mathieuj ( f1999[ > 
hypothesized that about 1 Myr ago, a supernova disrupted the 
parent molecular cloud and created the giant HII region seen 
to date, which is maintained by the remaining population of 
massive stars in the region. The expansion velocity of the H II 
region is traced by observations of the dense molecular shell 
surrounding the ionized gas, and is estimated at v exp =16.5 km 


s -1 expanding from the point (l,b) = 095.8.-12.11 centered a t 
vlsr = 3.8 km s 1 in velocity space (Lang & Mashederi|1998). 

The main ionizing star A Ori lies right at the (projected) 
heart of the bubble, which seems inconsistent with its proper 
mot ion (fi a cos 5 =-0.34+0.60 mas and fig = -2.94+0.33 mas, 
van Leeuwen 2007) with respect to the dense, molecular 


shell |MaddaIena & Morris||l9871 Lang & Masheder 19981. 

We calculate the LS R space v elocity of A Ori following the 
method described in Cox et al.| < 2012 >, using its proper mo¬ 


tion 


(van Leeuwen|20Q7 1 , radial velocity (v ia( j = 30.1 km 


Gontcharov||2006|l, while correcting for the solar motion us- 
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Fig. 5.— Three-color image of the A Ori bubble. The color codes are the 
same as Fig. [III. Overlayed are contours of IRAS 60 fim, with intensity lev¬ 
els: maximum - 90% - 80% - 70% - 60% - 50% - 40% - 30%. The solid 
line marks the past trajectory of the star A Ori during 4 Myr given its cur¬ 
rent space motion. Open diamond symbols mark intervals of 1 Myr. The 
dashed lines are the extremes of possible space motion considering the errors 
in proper motion and distance to the star. The black triangl e marks the centre 
of expansion of the A Ori bubble jLang & Masheder|1998) . 


ing parameters from Cogkunog lu et al. 
distance of 450 pc to the star ([Dolan & 


(j2orn>. 

Mat fuel 


We adopt a 


athieu|2001(l. The 


calculated space motion of A Ori is i’ ra d = 18.2 km s 1 with a 
position angle of 330 degrees (from north-to-east) inclined by 
55 degrees in the plane of the sky. The trajectory is overplot¬ 
ted in Fig. [5] and reveals that it is unlikely that its current mo¬ 
tion has been imparted at birth, since A Ori would have then 
originated at the boundary between the molecular and ionized 
gas (the expected age of A Ori is ~3 - 5 Myr; see Hernandez 
|et al.|2 010 and references therein). However, its motion and 
distance from the center of expansion (Fig. [5} is consistent 
with the SN hypothesis that would have launched /I Ori on its 
current trajectory about 1 Myr ago. 

We recognize the arc-shaped mid-infrared emission around 
A Ori, which is the tell-tale signature of a dust wave, created 
by the interaction of dust entrained in a champagne flow with 
radiation pressure from the sta r, which is a common phen om- 
ena inside interstellar bubbles ( |Ochsendorf et al.|2014a|b) . In¬ 
deed, the radio continuum maps presented by |Reich| (| 1978 1 
reveal an emission gradient within the ionized gas that is ex¬ 
pected for a champagne flow (Tenorio-Tagle 1979). The emis¬ 
sion gradient and the location of the dust wave reveals that 
the interior of the A Ori bubble is venting into the surrounding 
medium towards the south-east. However, the overall distri¬ 
bution of the ionized gas is still roughly spherical (Fig. HI; an 
extended emission region that would accompany a developed 
H II region champagne flow is not detected. This may indicate 
that the champagne phase has just commenced and the bub¬ 
ble is on the verge of breaking out. Th e imp lications of these 
observations shall be addressed in Sec. 14.31 

The WHAM observations in Fig. |4] reveal that the A Ori 
H II region lights up at positive vlsr (Fig. |4|, extending long- 
wards of v'lsr > 80 km s 1 . In addition, elevated emission 
surrounding the A Ori region can be seen at the other veloc¬ 


Cowie et al. 


ity extreme of WHAM, i.e., -80 km 
s 1 , possibly tracing the same component from 

who detected faint HV gas features towards A Ori. 
/HM of the Her line in the A Ori region is ~80 km s~ 1 , 
much broader than the thermal linewidth (~ 10 km s 1 ; Sec. 


3.3.1 1 . Thus, both the HV gas and Her suggests higher veloc¬ 


ities than that observed expansion velocity of the molecular 
shell (v eX p = 16.5 km s _1 ; Lang & Masheder 1998). The dis¬ 


crepancy between velocities observed for the dense molecular 
ring and the ionized gas could indicate a fast blister flow (see 
above) that can lead to large line of sight variations, or a re¬ 
cent explosive event that has accelerated the gas in the bubble 
interior, but has not yet coupled to the dense molecular ring. 

3.3.4. GS206-17+13 and the Orion nebula 

Besides the obvious large-scale shells that are associated 
with Barnard’s Loop and A Ori, we note that the Orion com¬ 
plex harbors shell structures on all detectable scales. Be¬ 
low, we will discuss the degree-sized GS206-17+13 shell 
(Ehlerova & Palous 2005) and the shell s urrounding the Orion 
Nebula ( |O’Dell|2001 j|Gudel et al.|2 008 ), but on even smaller 
scales there have been detections of many regions within the 
Orion A cloud that appear to have shell-like structures, pre¬ 
sumably driven by Fill regions from massive young stellar 
objects ( [Bally et al.p987| |Heyer et al.|T992) . Both GS206- 
17+13 and the Orion Nebula emit at positive vlsr (Fig. HI 
with large line-widths (FWHM ~70 km s" 1 ) that can both be 
related to ionized (blister) flows as described below. 

IQchsendorf & Tielens ( 12015) 1 discussed the GS206-17+13 
shell which is directly adjacent to the OMC B, has dimensions 
of 2° x 4° (14 x 28 pc), and is approximately centered on the 
cr Ori star cluster. The formation of the shell is likely caused 
by the stellar winds of the Orion OB lb subgroup that formed 
the protrusion out of the OMC B. The bright Hu emission 
that is projected inside part of the shell (the IC 434 emission 
nebula) is not related to the formation or history of GS206- 
17+13, as this gas originates from the champagne flow driven 
by the ionizing flux of the cr Ori cluster that has entered the 
pre-existing bubble cav ity and is now approaching OMC B 
(Ochsendorf & Tielens 12015 1. Ionization of the cloud edge 
sets up a pressure discontinuity between the cloud and inter 
cloud medium, and drives a shock front into the low-density 
bubble, while a rarefaction wave travels into the OMC B 
cloud. In panel 2a of Fig. [2] the shock front is visible at 
r s , whereas the rarefaction wave is located at r lw . Between 
the shock fronts, a champagne flow of ionized gas is set up 
(Tenorio-Tagle) 1979|), flowing from r rw towards r s . 

The density structure of the champagne flow of the IC 434 
region is well-ch aract erized by an exponential density gra¬ 
dient (Ochsendorf & Tielens 2015), reaching from tin = 35 
cirT 3 at the IF on the OMC B cloud surface, to 3.5 cirT 1 
at the location of cr Ori (r*; Fig. [2]) where the flow reaches 
~30 km s _1 (Mach number Al - 3 for a 10 4 K gas), which is 
the typical maximum velocity reached for champagne flows 
(Tenorio-Tagle 1979). Therefore, the ionized gas will not ac¬ 


celerate much further (also reflected by the flattening of the 
emission profile behind r*; Fig. [2]) and the density in the 
flow remains roughly similar up until the shock front, i.e., 
pi - 3.5 cm -3 . Then the pre-shock density po in which the 
shock front is moving can be determined through p\/M 2 = 
po (Bedijn & Tenorio-Tagl e|l98T I, from which we infer po = 
0.4 cm A which is far below tEat observed for electron den¬ 
sities within evolved Galactic Hu regions of similar size (> 
50 cirT 3 ; Paladini et al. 2012), however typical for densities 
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observed within the WIM ( ]Tielens|2005) l. Moreover, the ex¬ 
istence of the cometary clouds L1617 and L1622, which are 
located behind Barnard’s Loop as measured from the Orion 
OB association and have sharp I-fronts on their sides facing 
the Belt stars ( Bally et al.|2009j > reveal that the ionizing flux 
of the Belt stars must escape the GS206-17+13 shell, from 
which we conclude that the shell is incomplete and the bubble 
has likely ‘burst’. Hence, we posit that the bubble interior has 
been emptied in its past, and its contents was channeled into 
the Orion-Eridanus superbubble cavity. 

The Orion nebula is a prototypical blister Hn region ere 
ated by the young stars of the Trapezium (<1 Myr; Hillen 


brand|| 1997| [Bally 


2008) with 9 l Ori C as the dominant ion 
izing source. The ionized gas is streaming fro m the near-side 
of the OMC at about 10 km s 1 (O’Dell 2001) and accelerates 


into the surrounding medium. Recent velocity measurement 
of 9 l Ori C reveal that the star may be moving_away from the 

and 


O’Dell et al. 2009 


OMC at a velocity of 13 km s _1 (see 
references therein), which would imply a dynamicaltlmescale 
of the Orion nebula of only 1.5 x 10 4 yr, since the star is cur¬ 
rently located ~0.2 pc away from the IF. Indeed, the measured 
EM reveals that the thickness of the ionized gas layer is only 
0.13 pc ( |Wen & Q’Dell||T995 1 assuming a constant density, 
which is known to be ~6X 1 0 cm~ 3 at the IF ( O’ Del 1 1200 1 1 . 
X-ray emission from the ~10 6 K gas of the shocked stellar 
wind of 9 Ori C ( |Oudel et al.|2008 1 has been detected in the 
Extended Orion Nebula (EON) and confirms the leakage of 
the (hot) gas towards the south-west, with the Orion-Eridanus 
superbubble as the most likely outlet of the flow of material. 

4. DISCUSSION 


Cowie et al. (19791 interpreted neutral absorption lines at 
|vlsrI Si 20 km s ' towards stars in Orion as a thick, dense 
shell of swept-up material, partially ionized on the inside by 
the Orion OB association. This ionized component would 
be evident as the IV gas and Barnard’s Loop. However, the 
LV and IV gas are clearly separated in velocity space (Welty 
et al.] [2002| > towards f Ori at (/, b) = (206.5,-16.6). This ve¬ 
locity separation makes it unlikely that LV and IV trace the 
same component of the Orion-Eridanus superbubble. In this 
respect, Bally| ( |2008| > already noted that the velocity disper¬ 
sion of the ensemble of atomic and molecular clouds in Orion 
ranges over > 20 km s -1 . Therefore, we attribute the low- 
velocity neutral gas to clouds in the Orion region instead of 
a dense, neutral shell surrounding Orion. In addition, we 
have revealed the absenc e of a dense neutral shell surround¬ 
ing Barnard’s Loop (Sec. |3 ,2[ >, and have shown that the Loop 
is part of a (separate) complete shell structure that is spatially 


correlated with HV gas detected towards Orion (Cowie et al. 

a 


1979] Welty et al.|2002]> as discussed in Sec 
The above described results call for a thorough revisit of 
the structure, history, and subsequent evolut ion of the Orion- 
Eridanus superbubble. Below in Sec. |4.1.1| we will connect 
the different gas tracers to the components of the superbubble, 
and combine this with our observational work in order to de¬ 
rive a n ew and impro ved pi cture of the Orion-Eridanus region 
in Sec. |4.1.2|& Sec. |4.1.3| We relate our findings to the stel¬ 
lar content of Orion OB 1 in Sec. 4.2 and present a scenario 
for t he fu ture evolution of the Orion-Eridanus superbubble in 
Sec. 1431 


4.1. Connecting the dots: an updated picture of the 
Orion-Eridanus superbubble 


4.1.1. Probing the superbubble through absorption lines 
Intermediate-velocity gas components are detected over a 

J a SO Cm o “ 1 ^ ,,_^ OO Cm o ~ 1 


(Cowie & York 1978 

Cowie et al. 1979, Huang et al. 1995, 

Welty et al. 2002). I 

lowever, with high spectral resolution, 


et al.](j2002f were able to locate the bulk of the IV gas towards 

Ori in two single components at vi sr = -35.3 km s 1 and vi sr 
= -31.6 km s _I . 

We link the IV gas to the outer superbubble wall through 
several observations. First, the bulk of the IV gas lies within 
the velocity space reported for the superbubble wall in Hi (- 
40 km s _I < vlsr 55 40 km s -1 ; Brown et al.|l995 1. We note 
that by connecting the (ionized) IV to the superbubble wall, 
we assume that the outer wall towards ( Ori is predominantly 
ionized along the line of sight. This might connect with our 
result that the eastern part of the superbubble wall does not 
exhibit a neutral shell, as was observed in reg ion 1 of Fig. 
|T| and will be further discussed in Sec. |4.1.3 Second, the 
total Hll column density (log[N(Hn )] = 18.44 cnL 2 ) and 
the post-shock conditions of the IV gas towards ( Ori (n \ = 
0.16 cm 3 , T\ = 9e3 K; Welty et al.|2002 1 imply an ionized 
shell thickness of ~5.6 pc, in close agreement with the mea¬ 
sured thickness of 1 degree or ~1 pc of the Her filaments along 
regions 1 - 5 in Fig. IT] Third, the IV gas is widespread, ex¬ 
tending well outside Barnard’s Loop and encompassing the A 
Ori region, tracing a distinc t component over an area larger 
than the HV gas (Sec. |3.3.1| l. Even though the definite extent 
of the IV gas should be established through an absorption line 
study that include sightlines covering the entire region shown 
in Fig. [4] here we argue that that the IV gas traces a structure 
significantly larger than the dimensions of the Orion-Eridanus 
superbubble described in previous works (|Reynolds & Ogden] 


[1979][Heiles et al.|l999irBallyi2008llPon et al.l2014af~ 

The above results reveal that the IV gas may be linked to 
the outer superbubble wall, whose n eutra l component was de¬ 
fined in Brown et al. (1995b. In Sec. |3.3| we have argued that 
the HV gas traces the fast expansion of the Barnard’s Loop 
bubble through the correlation between the sightlines that ex¬ 
hibit HV features and the morphology of the Barnard’s Loop 
bubble from Ha (Fig. [4]). A similar result was obtained for 
the A Ori bubble. Both the Barnard’s Loop bubble and the A 
Ori bubble may be connected to recent SN explosions (Secs. 


3.3.2 & 3.3.31 that could be at the origin of the HV gas. Fi¬ 
nally) the LV gas is attributed to the ensemble of clouds in 
Orion ( |Bally||2008| ). As a caveat, we note that these conclu¬ 
sions are based on the study of a limited amount of sightlines 
that may not fu lly reflect the complexity of the Orion region 
(see Sec. |4. \A\. Nonetheless, based on the currently available 
data, we hypothesize that the different gas tracers trace sepa¬ 
rate components of the superbubble, summarized in Table [2] 

4.1.2. Extent of the outer superbubble wall 

The exte nt and connection of the IV gas to the superbubble 
wall (Sec. |4XI) imply that the size of the Orion-Eridanus 
superbubble may be larger than previously thought. Indeed, 
Barnard’s Loop is part of a separate closed bubble structure, 
suggesting that it is not associated with the large-scale ex- 
panding superbubble structure defined in earlier works (e.g., 


Reynol ds & Og den|1979 Heiles et al.|1999| Bally|2008] Pon| 

et al.||2014a| l. Here, we argue that the outer Her filament to¬ 

wards tEe"east in Fig. [Tjoffers a prime candidate for represent¬ 
ing the outer superbubble wall. 
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TABLE 2 

The connection between observed gas velocities and components within the Orion-Eridanus region 


Component 

vlsr 
(km s -1 ) 

no 

(cm -3 ) 

n\ 

(cm -3 ) 

log[(V H ] 

(cm^ 2 ) 

ion. state 

Ref. 

Tracing 

Low-velocity (LV) 

-20 to 20 

- 

- 

20.43 

neutral 

1 

Atomic/Molecular clouds in Orion 

Intermediate-velocity (IV) 

-80 to -20 

- 

0.16 

18.44 

ionized 

1,2, 3,4 

Orion-Eridanus superbubble wall 

High-velocity (HV) 

-120 to -80 

3 x 10^ 3 

0.1 -0.2 

17.88 

ionized 

1,2 

Nested shells (Barnard’s Loop bubble, A Ori) 


Note. — Listed quantities are (if known): observed velocity extents, vlsr; pre-shock densities, no; post-shock densities, n i ; total hydrogen column densities, 
Nn ; ionization state, and the components of the Orion-Eridanus superbubble that is traced by the specifi c ga s velocities. The (colu mn ) densities quoted ar e 
observed in the line of sight towards f Ori. References are (l): |Welty et al.|f2002) , (2): |Cowie et al.|jl979| . (3): |Cowie & York|(l 978) (4): |Huang et ai.|(i~995f 


We have concluded that Barnard’s Loop is completely ion¬ 
ized by Orion OBI and that ionizing photon s a re abl e to 
penetrate through. Moreover, in Secs. |3.2.1| and |3.2.2[ we 
have shown that Orion OBI can provide the necessary ion¬ 
izing photons to illuminate the Her emission seen in Fig. [T] 
indicating that the ionization front lies well outside the pre¬ 
viously defined bounds of the Orion-Eridanus superbubble, 
which may be possible for superbubble walls in general (Basu 


et al. |1999) >. However, the large-scale, coherent Her filament to 
the north-east is indicative of a limb-brightened shell, rather 
than a mere haphazard ISM structure which happened to be 
illuminated by ionizing photons piercing through Barnard’s 
Loop. 

We note that often, the extent of superbubbles is traced 
through diffuse X-ray emission. The photo-electr ic cross sec- 
tion for 0.25 keV photons is ~4 x 10 --1 cm 2 H _1 (Morrison & 


McCammon|1983|l. By integrating the total foreground emis¬ 
sion from the LAB survey (Kalberla et al.||2005 


i.e., -400 

< Vlsr < 0 km s -1 (which includes the near-side of the su 


perbubble shell at velocities -40 < v'lsr < 0 km s 1 ; [Brown 
jet al,|1995) , we infer that the 0.25 keV emission is confined to 
foreground Hi column densities of (V H S3x 10 19 cm ~ 2 (the 
superbubble cavity), where the optical depth for 0.25 keV X- 
ray photons is r x <0.1, while emission at 0.25 keV is absent 
from regions of Nh > 2 x 10 20 cm ~ 2 (i.e, Barnard’s Loop, and 
the east- and westside of the superbubble near the Galactic 
plane), where r x > 0.6. Thus, the clear anti-correlation of the 
0.25 (and 0.75) keV emission (Fig. [T]) with foreground HI lim¬ 
its the use of this tracer for many parts of the Orion-Eridanus 
region, especially towards the Galactic plane. 

The association of Barnard’s Loop with a closed separate 
bubble, the IV gas that traces a region larger than the pre¬ 
viously defined bounds of the Orion-Eridanus superbubble, 
and the mere presence of the large-scale, coherent Her fila¬ 
ment to the north-east, has led us to believe that regions 1-5 
(Fig. trace the actual outer shell of the Orion-Eridanus su¬ 
perbubble. Projected in the interior of the outer shell reside 
the Barnard’s Loop bubble, the A Ori bubbl e, and t he smaller- 
scale expanding bubbles discussed in Sec. |3.3.4| This indi¬ 
cates that the superbubble in its entirety consists of separate 
structures, or successive nested shells , possibly connected to 
a series of SN explosions originating from Orion OBI over 
the past 15 Myr ()Bally|2008( Secs. |4“L4|&|4T2]l. While hot 
plasma could exist between the Ha filament and Barnard’s 
loop, the detection of the accompanying X-ray emission is 
hampered because of foreground absorption. Our hypothesis 
sprou ted from s tudyin g the energetics of gas and dust (Sec. 
|3.2.1| and Sec. |3.2.2|) and their morphological appearances 
(rtg- 0 - combined with the connections between the differ¬ 
ent gas tracers with the components of the superbubble (Table 
[2|. In the next subsection, we will discuss the temperature of 
the superbubble interior gas, as well as the global structure of 
the outer wall in order to substantiate our claim regarding the 


structure and extent of the Orion-Eridanus superbubble. 

4.1.3. Temperature of the superbubble interior and structure of the 

outer wall 


The evidence discussed in Sec. I4.1.2l indicate that th e oute r 
superbubble wall lies outside of Barnard’s Loop (Sec. |4. 1.2|), 
possibly extending towards the outer Her filament to the north¬ 
east (Fig. HI. Here, we quantify the temperature structure of 
the superbubble interior. First, we derive that the Barnard’s 
Loop bubble is moving supersonically towards the east, from 
which we conclude that the temperature in this direction can 
not be as high as >10 6 K. After comparison with the cooling 
timescale for a gas at <10 6 K, we infer that the medium be¬ 
tween Barnard’s Loop and the outer Her filament must have 
cooled to ~10 4 K, appropriate for a gas photo-ionized by 
Orion OBI. We argue that mass-loading may cause tempera¬ 
ture gradients within superbubble interiors, place our findings 
in the context of superbubble evolution, and discuss the struc¬ 
ture of the outer wall. 

The temperature structure of the superbubble interior can 
be estimated as follows. First, we define the isothermal sound 
speed, c s = ( kT/panw )^ 2 , where p is the mean mass per hy¬ 
drogen nucleus (= 0.61 for a fully ionized medium). Towards 
the west of the Barnard’s Loop bubble, the temperature of the 
superbubble interior traced by the X-ray emitting gas is esti¬ 
mated at T x = 2.1 x 10 6 K (Burrows et al. 1993 i, giving c s = 
170 km s 1 . We compare this w ith the Barnard’s Loop bubble 
expanding at -100 km s 1 (Sec. 3.3.2 1 , which implies that the 
remnant of the recent SN explosion (Secs. IS & |43T) is 
well within the radiative expansion phase (see below). During 
this phase, the velocity of the SNR decelerates with time as 
cc r 5n (Tielens 2005) if we assume that the ambient density 
in which the SNR propagates is distributed homogeneously. 
Adopting this time-velocity relation, the adiabatic expansion 
phase (v exp > 250 km s 1 ) ended just below 10 5 yr, and the 
bright crescent of Barnard’s Loop has been moving subsoni- 
cally for ~0.2 Myr (in case c s = 170 km s -1 ). When moving 
subsonically, sound waves will travel ahead of the shell and 
the structure would dissolve and merge with the ISM within 


a fade-away timescale of (McKee & Ostriker 1977 

|20TT| l 


Draine 


Tf ade = 1.87 E\ 


0.32 -0.37 I 


51 


n 


t c s \-V 5 

(=) Myr - (5) 


where £gi is the kinetic energy in numbers of 10 M erg, and 
n the ambient hydrogen number density. The current kinetic 
energy of the eastern half of the Barnard’s Loop bubble mea- 


sures £ 5 ! ~ 0.34 (Sec. 4.3.1 1 , part of which may already have 

been thermalized (Mac Low & McCray 

1988). However, the 

total kinetic energy delivered by a SN ex 

fiosion typically does 
), or £51 ~ 0.5 for the 

not exceed £51 ~ 1 (Veilleux et al. 2005 


eastern half of the bubble. With £51 = 0.5, and parameters for 
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Fig. 6. — Schematic of the Orion-Eridanus superbubble and several of its 
major components. See Fig. ITlfor the nomenclatures of the shown structures. 
The proposed outer shell of tne Orion-Eridanus superbubble (dashed circle), 
blown by a series of SNe from an old subpopulation of the Orion OB asso¬ 
ciation, is traced by an intermediate-velocity (IV) shock that sweeps up the 
ambient ISM (gray area). Towards the west, the superbubble is surrounded 
by a shell of neutral (solid black line) swept-up material, whereas this shell 
is completely ionized (dashed gray line) towards the east. The temperature 
in the superbubble interior seems to be lower towards the east (~10 4 K) in 
the direction of the Galactic plane, compared to the X-ray emitting gas in the 
west (>10 6 K). The outer shell encompasses several ‘nested’ smaller shells 
or bubbles (solid gray circles, denoted in boldface) that are more recently 
formed, such as the Barnard’s loop bubble, traced by high-velocity (HV) gas. 
Other examples of bubbles that may have been triggered by ongoing activity 
in Orion OBI are the A Ori bubble, GS206-17+13, and the shell surrounding 
the Orion Nebula cluster (ONC). 


the X-ray emitting gas of the superbubble interior towards the 
west (T x = 2.1 x 10 6 K, n x = 6 x 10 -3 ctrT 3 ; Burrows et al. 


1993 


Sec. 4.3.11, we calculate Tf a d e ~ 0.2 Myr, similar to the 
time Barnard’s Loop has been moving subsonically in case it 
is moving within hot gas at a temperature exceeding 10 6 K. 

The above derived fadeaway timescale reveals that if 
Barnard’s Loop would be moving to the east in an X-ray emit¬ 
ting gas such as that traced by ROSAT in the west part of the 
superbubble (Fig. [lj, it must be long in the process of merg¬ 
ing and dissolving with the hot gas. This actual process may 
cause the faintness of the Her e missio n of the Barnard’s Loop 
bubble towards the west (Sec. |3.3.2[ ), where the superbubble 
interior temperature _is known to be at a temperature of T x = 
2.1 x 10 6 K (Burrows et al. 1993; Fig. Jl}. Clearly, this does 
not affect the eastern part of the Barnarcrs Loop bubble in the 
same way, as exemplified by the bright crescent of Barnard’s 
Loop. Alternatively, there may be a density gradient within 
the superbubble interior that may cause the brightness differ¬ 
ence of the Barnard’s Loop bubble between its east and west 
side. Either way, the conclusion is that the Barnard’s Loop 
bubble moves supersonically towards the east, in line with the 
findings of Cowie et al.|(|l979jl, who argued that the HV gas 
(ass ociated with the near-side of the Barnard’s Loop bubble; 
Sec. 3.3.2[ ) traces a radiative shock. Consequently, the interior 
superbubble temperature towards the east can not be at >10 6 
K, as is observed for the X-ray emitting gas towards the west. 

Temperature differences within a superbubble interior can 


originate from thermal conduction and evaporation of the 
swept-up superbubble shell. This mechanism may cool the 
interior of the superbubble by increasing the interior density 
( [Weaver et al.|[T977] |Mac Low & McCray||1988] l such that 
radiative losses become important. Champagne flows and 
photo-ablation of molecular clouds enclosed in the superbub¬ 
ble can inject mass into the hot plasma, further lowering the 
temperature compared to that predicted by idealized models 
(a full discu ssion on the mass-loading mechanisms is deferred 
4.3 i. Gas below temperatures of 10 6 K cools rapidly. 

10 -3 


to Sec. 


-30 


At an interior density of n x = 6 x 10 3 cm 1 (Burrows et al. 
1993), the cooling timescale for a gas at 10 5 K equals r coo i ~ 


-50 


kT/nA(T) ~ 10 s yr, where A(T) ~ 5 x 1(L 22 erg cm J s _1 is the 
average value of th e cooling function for gas at temperatures 
10 4 K < T < 10 6 K (Dalgarno & McCray 19721. Mass-loading 
of the interior will shorte n the coo ling timescale. The detec¬ 
tion of N V gas (Welty et al. 2002) towards f Ori traces ~10 5 
K gas and confir ms the presence of a thermally conductive 
layer ( |Chu 2008 1 , providing information on the temperature 
structure of the superbubble wall and the interior temperature, 
but it is unclear if this conductive layer traces the outer super¬ 
bubble wall, or the interior shell associated with the Barnard’s 
Loop bubble. Furthermore, it is uncertain if this layer extends 
towards the outer Ha filament given the lack of coverage of 
absorption line studies in this region (Fig. [4}. In any case, the 
cooling timescale is short compared to that of the estimated 
age of the superbubble (< 8 Myr; Brown et al. 11994 ) and the 
supernova rate in Orion OBI (1 - 1.5 Myr -1 ; Bally 2008] Sec. 
4.3.2| i that may re-heat the interior gas (see below and Sec. 
4.3[ ). Consequently, if the temperature drops below 10 6 K, the 
medium between Barnard’s Loop and the outer filament will 
settle at a temperature of ~ 10 4 K, as it is kept photo-ionized 
by O rion O BI by photons piercing through Barnard’s Loop 
(Sec. |3.2.2|). Depending on the internal density, the medium 
between Barnard’s Loop and the outer filament may not be in 
pressure equilibrium with the (higher-density) outer Ha fila¬ 
ment. 

In the early stages of its evolution, an outer superbubble 
wall is defined by a shell of swept-up material, the expan¬ 
sion of which is driven by an overpressurized interior medium 
(Tielens |2005Jl. The superbubble initially expands through 
the Sedov-Taylor solution, but as the interior medium ex¬ 
pands adiabatically, its pressure will decrease and the expan¬ 
sion slows down. The radiative phase commences when the 
expansion velocity decreases_to ~250 km s or a post-shock 
temperature of ~10 6 K (Tielens 2005). At this point, radia¬ 
tive losses will become important, and the interior tempera¬ 
ture (and pressure) will further drop. When the interior pres¬ 
sure reaches that of the ambient medium, the expansion of 
the superbubble wall will coast at constant radial momentum, 
thus slowing down as it sweeps up material from its surround¬ 
ings. When the expansion velocity reaches the local sound 
speed, the superbubble wall will disappear and merge with 
the ISM within a fadeaway timescale (Eq. |5j. The interior 
pressure may even drop below that of the ambient medium 
if the cooling timescale (~10 5 yr) is shorter than the pres¬ 
sure timescale or sound-crossing time. Currently, the dis¬ 
tance between Barnard’s Loop and the outer Halpha filament 
is roughly ~100 pc, such that the pressure timescale can be 
anywhere up to ~10 Myr for a 10 4 K medium, depending 
on the exact moment the region cooled towards this tempera¬ 
ture. However, the medium between Barnard’s Loop and the 
outer filament must have cooled within the time between sue- 
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cessive expanding shells (regulated by the SN rate; 1-1.5 
Myr -1 ; Bally 2008). Thus, it may well be that the pressure in¬ 


side the interior cannot adjust rapidly enough to the change in 
temperature and, in case no mass-loading occurs, it can reach 
a lower pressure compared to the ambient galactic medium. 
In case of an underpressurized interior, the expansion of the 
superbubble wall may halt and even reverse its movement if 
the pressure difference between the interior and the ambient 
medium is strong enough. 

Given that the medium between Barnard’s Loop and the 
outer Her filament is at 10 4 K and we are not able to constrain 
its density, one may question if this medium in fact represents 
the WIM or WNM. In this case, the outer Hu filament is not 
the superbubble wall, but merely an ISM feature illuminated 
by photons escaping from Barnard’s Loop. However, we have 
revealed t hat Ba rnard’s Loop is part of a separate closed struc¬ 
ture (Sec. |3.3.2|), and argued that the superbubble wall extends 
beyond Barnard’s Loop (addressed in Secs. |4. 1.1 1 & |4.1.2|i. 
Surely, if one accepts that Barnard’s Loop is expandi ng at 10 0 
kmsj 1 and represents a SNR of age 0.3 Myr (Secs. 


3.3.2 


& 


4.3.1 1), it may not be surprising that the superbubble wall lies 


further outwards as measured from Orion OBI, given that star 
formation has occurred over 15 Myr in Orion and that 10-20 
SNe have gone of in its past (jBally||2008). Here, we argue 
that successive SN explosions may lead to the appearance of 
successive shells. In this scenario, the large-scale, coherent 
Hu filament to the north-east represents a prime candidate to 
trace the actual outer superbubble wall. We defer a further 
discussion on the idea of successive shells and the evolution 
of the outer superbubble wall to Sec. |4.3.2| 

We note that the outer Hu filament (Fig. |T|) does not con¬ 
tain a neutral shell that is expected for a sweeping superbubble 
wall (Fig. [2]&[3]). The lack of a neutral shell may simply sug¬ 
gest that the swept-up ISM did not contain enough material 
to fully absorb the ionizing photons from Orion OB 1. Alter¬ 
natively, the east side of the superbubble, which is racing up 
the density ramp towards the Galactic plane, may not be mov¬ 
ing supersonically anymore and has changed character from 
a shock wave to a sound wave, rendering it unable to sweep 
up additional mass. Using Eq. |5]a nd E $i ~ 10 (the total ki- 
netic energy exterted Orion OBI; Brown et al.j |1995j [Bally 
2008]), we calculate for the east side of the bubble expand¬ 
ing towards the WNM of the Galactic plane (n ~ 0.5 cm~ 3 , 
c s ~ 10 km s -1 ), Tf a de — 4 Myr, well below the estimated age 
of the superbubble (< 8 Myr; Brown et al.|| 19941) . Thus, we 
may be observing the east side of the superbubble as it is in 
the process of dissolving, already having lost its neutral shell 
and currently being completely ionized before it merges with 
the local material of the ISM. Eq. [5] does not apply to the 
west side of the superbubble as its expansion velocity (traced 
through Hi ; ~40 km s~ ! ; Brown et al. 1995) implies that it 
moves supersonically through the WNM/WIM. However, for 
the ionized east side of the superbubble the expansion veloc¬ 
ity is not constrained, and the question whether or not it is 
moving supersonically remains unanswered. 

4.1.4. Summary: the structure of the Orion-Eridanus superbubble 

We now have obtained the complete set of tools to derive 
an updated picture of the Orion-Eridanus region which we 
show in Fig. [6] Our results demonstrate that the superbub¬ 
ble consists ofseparate structures or nested shells, superim¬ 
posed along line of sight. In particular, we h ave recognized 
a separate bubble (amongst others; Sec. |3.3| ), evident as the 
HV gas and the Barnard’s Loop bubble in Ha (Fig [4]t. This 


bubble expands inside the pre-existing cavity of the Orion- 
Eridanus superbubble, given its pre-shock density, its embed¬ 
ded nature in the Orion-Eridanus region (Fig. [6]), and its likely 
origin through feedback from Orion OBI. The realization 
that B arnard’s Loop is not part of the superbubble wall (Sec. 
3.3.2), the IV gas that seems to trace a region larger than the 


previousl y defi ned bounds of the Orion-Eridanus superbub¬ 
ble (Sec. 4.1.1 1 , and the presence of the large-scale, coherent 


Her filament to the north-east, has led us to believe that the 
extent of the Orion-Eridanus superbubble is larger than pre¬ 
viously throught. We have argued that the temperature of the 
superbubble interior between Barnard’s Loop and the outer 
Her filament to the north-east can not be as high as that ob¬ 
served for the X-ray emitting gas towards the west (Burrows 
et al. 1993) and, therefore, must have settled at l(r K (Sec. 


4.1. 3). The expansion of this outer Her filament may not be 


pressure-driven anymore. While this is an inevitable outcome 
during late stages of superbubbles and SNRs (Tielen s|2005 ), 
the interior pressure drop may be accelerated through t he ef - 
fects of mass loading and radiative energy losses (Sec. |4.3b. 
The observations (Fig. IT]) indicate that the outer superbubble 
wall is surrounded by a shell of HI towards the wes t, wher eas 
it may be completely ionized towards the east (Sec. |4. 1.3). 

We stress that the association of gas velocities with com- 
ponents in the superbubble (Table |2]) are largely based on the 
analysis of a dozen of sightlines throughout the Orion side 
of the superbubble (Fig. |4]. It is not cle ar how much mass 
is associated with these velocities (Sec. |4.3.1| >. Moreover, 
some of the gas tracers are not seen in particular lines of sight. 
This may be because we are tracing foreground stars, such as 
the sightline towards y Ori located at 77 pc (van Leeuwen 


2007), which could provide a constraint to the elo ngation of 


the bubble towards the Local Bubble (e.g.. Bally 


et al. 


2008 Pon 


|2014a), or incomplete spectral coverage (|Cowie et al. 


1979 . Alternati vely, the d istribution of the gas components 
might be patchy (Cowie e t al,|1979) . Surely, the connections 
laid out in Table[2]may be oversimplified, as the Orion region 
is highly complex, containing dozens of filamentary clouds, 
and must bear the marks of multiple SNRs ( |Bally|2008] ) that 
are breaking out of a region that contains a lot of substruc¬ 
ture. Further observations are needed to firmly constrain the 
extent of the distinct gas components in order to derive the 
detailed structure of the Orion-Eridanus superbubble, and to 
recognize the existence of conductive layers that would pro¬ 
vide more insight in the temperature distribution throughout 
the superbubble interior. 

The observations of the gas and dust have revealed that 
Orion OBI can be the source of ionization of the Ha fila¬ 
ments encompassed within the entire region shown in Fig. [T] 
(Sec. |3.2) . We have connected the IV gas with the outer su¬ 
perbubble wall. We note that the pre-shock density of the IV 
is not known. However, given a post-shock density of 0.16 
cm~ 3 (Welty et al. 2002) and shock velocity of 40 km s 1 
( Brown et al.||T995) , a shock moving through the WNM (T 
- 8000 K, c s ~ 10 km -1 ) would have a pre-shock density of 
order -0.01 cm~ 3 (J-type shocks increase densities by AT 2 ; 
|Tielens||2005) , which is somewhat low for the WNM in the 
Galactic plane. Higher densities for the Orion-Eridanus shell 
of 1 - 5 cm~ 3 are quoted towards the E ridanus filaments in the 
west (e.g., jReynolds & Ogden]|1979) , indicating that the ex¬ 
pansion occurs somewhat anisotropic ally, which can also be 
inferred by appreciating the complex morphology of the su¬ 
perbubble. In any case, the shock velocity of the IV gas (~ 
40 km s 1 ) is not fast enough to ionize hydrogen: therefore. 
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the Eridanus filaments (and the IV gas) are photo-ionized by 
Orion OB 1, which was alread y noted by [Reynolds & Ogden 
(19791 and Welty et al. (2002 1 . We note that Pon etal. ”. 


explored other mechanisms of ionization, but concluded that 
Orion OBI most likely causes the observed ionization state 
of the IV gas. When the spherical geometry from Fig. [6] is 
adopted for the outer superbubble wall, and assuming a di¬ 
ameter of the supershell of 45 degrees (equivalent to a radius 
of 160 pc) with expansion velocity of ~ 40 km s -1 , the dy¬ 
namical age of the shell is ~ 4 Myr, consistent with the age 
of the Orion OB la subgroup, minus a main sequence lifetime 
of 3 - 5 Myr for its most massive stars ( Bally||2008 ). Natu¬ 
rally, stellar winds during the first few Myr of the lifetime of 
the OB association may have contributed to the formation and 
expansion of the superbubble wall as well. 

In conclusion, considering that star formation had been oc¬ 
curring for 10-15 Myr in Orion OBI and th e presence o f 
at least four distinct subgroups ( jBlaauw |1964| |B ally |[2008 1 , 
it may not be surprising that the Orion-Eridanus superbub- 
ble consists of a set of nested shells. Our results connect 
well with recent work from |Schlafly et ak ( 2015] ), who show 
that the Orion molecular clouds are in turn part of an ancient 
ring of dust. Possibly, the ring represents the remainder of a 
large bubble that pushed the Orion molecular clouds out of the 
Galactic plane towards their current location. Still, the hierar¬ 
chy of bubbles does not end at the possible progenitor of the 
Orion-Eridanus superbubble defined by |Schlafly et al.|(2015[ i. 
Indeed, it has been argued that the Gould’s Belt that include 
all nearby OB associations are associated with a large expand¬ 
ing ancient (30 - 60 Myr) supershell known as Lindblad’s ring 
(see|Bally|2008| and references therein). The Gould’s Belt of 
stars and OB associations may therefore represent sequential 
star formation when the Lindblad’s ring cooled and collapsed 
to form dark clouds out of which the evolution of the Orion 
region commenced. 

4.2. The nested shells and their relation to Orion OBI 

The recent discovery that the Orion molecular clouds may 
be part of an ancient ring of dust (Schlafly et al. 20151, the un- 
co vering of an extended foreground population towards Orion 
A ( Alves & Bouy|20r2] |Bouy et al.'1 2014) , and the mysterious 
origin of some of the"brightest stars in Orion (a Ori, J3 Ori, 
and k Ori; Bally 2008) illustrates that our current understand¬ 
ing of the Orion region is still limited. |Bally|(]2008J) proposed 
that the bright stars may be part of a foreground group about 
150 pc in front of Orion OBI, and together constitute the three 
most massive members of this population still present to date. 
Alternatively, some of the massive stars from Orion OBI may 
have been ejected as runaways (space velocities in excess of 
> 40 km s 1 ) in the direction of the Local Bubble before ex¬ 
ploding as SNe: up to 2 5% o f the OB stars end up as run¬ 
aways (Zinnecker & York el|2007) . Such stars could be linked 
to the odd g eometry of the Orion-Eridanus superbubble ( |Pon| 
|et al.|2014aj), which appears to be elongated towards the Local 
Bubble. 

One must keep in mind that most of the distance determina¬ 
tions towards the Orion stellar populations are based on pho¬ 
tometric indicators (e.g. Brown et akj| 1994) . Although the 
Hipparcos mission determined parallaxes for many stars in 
the Orion region, the distance to the young populations (~ 
300 - 400 pc) translates to Hipparcos parallax errors of 30 - 
40% (given a precision of 1 milliarcsecond; Perryman et al. 
1997). Thus the detailed distribution of the stellar popula” 


(see also de Zeeuw et al.|l999 1 . In this re spec t, the Gaia mis¬ 
sion ( de Bruijne|| 2012 [ Berryman e t al.|[ 200 l i is expected to 
dramatically improve the distance and age determinations of 
the stars and subgroups in Orion. The entire Orion region 
lies well within the range where Gaia parallaxes will be ac¬ 
curate to 10 % or better over a large range of spectral types. 
This will allow to constrain the relationship between the stel¬ 
lar content and the nested shells seen in the Orion-Eridanus 
superbubble, and to address the importance of triggered star 
formation ( [Elmegreen & Lada|1977] l in Orion, for which one 
needs precise measurements of proper motion, radial veloci¬ 
ties, and ages. 

4.3. Evolution of the Orion-Eridanus superbubble 

In this section, we will exploit the results drawn from the 
previous sections to derive an updated picture of the evolu¬ 
tion of the Orion-Eridanus superbubble, regulated by a bal¬ 
ance between processes that mass-load the superbubble inte- 
rior, contro l ling its temperature an d density structure (Cowie 
et al j |1981 j Hartquist et al.1 1986) , and feedback from Orion 
OB 1 that sweep up and collect the injected inter ior mass in 
expanding shells before thermalizing or cooling (Mac Low &| 
McCr ayfl 988[ >. These processes determine the momentum of 
the superbubble wall and the thermal pressure of the interior, 
which drives the superbubble expansion out into the ambient 
Galactic medium. 


First, in Sec. 4.3.1 we summarize masses and energetics 
as observed from different components of the Orion-Eridanus 
superbubble. From this, we conclude that the Barnard’s Loop 
bubble is a SNR (asso ciated with a SN explosion that occurred 
3 x 1 0 5 yr ago; Sec. 3.3.1 Cowie et al. 1979), and argue in 
Sec. |4.3.2| that the SNR has swept through the inner parts 
of the superbubble to collect mass that has been introduced 
through champagne flows and thermally evaporating clouds. 
Behind the Barnard’s Loop bubble, champagne flows from IC 
434, the Orion Nebula, and A Ori continue to mass-load the 
superbubble interior until a next SN explosion will start the 
cycle all over again. The net effect is a gradual disruption of 
the star-forming reservoir and the subsequent outward trans¬ 
portation of processed material, where it will ultimately pres¬ 
surize the superbubble interior, or add momentum to the outer 
superbubble wall. We will argue that at the current rate, the 
OMCs will be able to power the expansion and evolution of 
the Orion-Eridanus superbubble for another 20 - 30 Myrs, be¬ 
fore they are depleted and run out of steam, at which point the 
superbubble will disappear and merge with the surrounding 
ISM. 

4.3.1. Masses and energetics of the Orion-Eridanus superbubble 
Elot gas (> 10 6 K; Table 3} from the superbubble interior 


has been characterized by 
ferred a temperature of T x 


__!P,_ 

Burrows et al.| (1993), who in- 

2 x 10 6 K and a total luminos¬ 
ity at X-ray wavelengths of L x » 8.0 x 10 35 erg s 1 , when 
assumed that the X-ray emitting region is spherical and its 
actual ‘hook-shaped’ appearance results after foreground ab¬ 
sorption by MBM 18 and MBM 20 (Fig. [T[ |Burrows et al.| 
1993). The authors deduce the following parameters for the 
hot gas: when 400 pc is taken as the distance towards the su¬ 
perbubble, the radius of the X-ray emitting region is r x = 70 


pc, its density is n x = 6 X 10 3 cm 3 (Burrows et al. 1993), and 


tion along the line of sight can, as of yet, not be determined 


contains a total mass M x = 230 M Q after adopting a spherical 
geometry. However, if we adopt the (larger) geometry from 
Fig. § and account for possible missing X-ray e mission be¬ 
cause of foreground absorption (Fig. [I] and Sec. 4.1.2 1 , the 
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TABLE 3 

Masses and energetics of the Orion-Eridanus superbubble 


M 

^exp 

Eth 

Gdn 

Form 

Reference 

(10 3 Mq) 

(km s -1 ) 

(10 49 erg) 

(I0 49 erg) 




Hot gas (ZIP 6 K) 


Superbubble interior 

1.8 

- 

62 

- 

Thermal energy 

Burrows et al. 

1993f 

Barnard’s Loop bubble interior 

(0.02) 


(0.65) 

- 

Thermal energy 

This work 


Ionized gas (10* K) 


Superbubble 

84 

15 (l > 

14 

19 

Swept-up shell 


Reynolds & Ogden 

1979 

Barnard’s Loop bubble 

6.7 

100 

1 

67 

Swept-up shell 

This work 






A Ori 

2-6 

30 

0.5 

3 

Champagne flow 


Reich 

1978}; 

van Buren 

1986) 

GS206-17+13 

3 

30 

0.5 

3 

Champagne flow 

This work 






IC 434 

0.1 

30 

0.02 

0.1 

Champagne flow 


Ochsendorf 

fc Tielens 

2015 

Orion Nebula 

0.02 

30 

0.003 

0.02 

Champagne flow 


Wilson et al. 

1997) 


Neutral gas (10 2 K) 


Superbubble 

250 

40 

- 

370 

Swept-up shell 

Brown et al. 

(1995 




GS206-17+13 

3.4 

8 

- 

0.2® 

Swept-up shell 

Ochsendorf 

& Tielens 

2015 


Orion Nebula (Veil) 

2.3 

2 

- 

0.01 

Swept-up shell 

van der Werf et al. 

20 

13); T 

iis work 


Molecular gas (10 K) 


OMC A 

105 

- 

- 

- 

Molecular cloud 

Wilson et al. 

2005 


OMCB 

82 

- 

- 

- 

Molecular cloud 

Wilson et al. 

2005 


/I Ori 

11 

16.5 

- 

6 

Swept-up shell 

Lang & Mas 

ieder 

1998 


Stars 


Ori la 

1.6 

- 

- 

- 

Stellar cluster 


Brown et al. 

1994 



Ori lb 

1.3 

- 

- 

- 

Stellar cluster 


Brown et al. 

1994 



Ori lc 

1.8 

- 

- 

- 

Stellar cluster 


Brown et al. 

1994 



Ori ld(ONC) 

1.8 

- 

- 

- 

Stellar cluster 


Hillenbrand 

1997 



Ori Id (NGC2024) 

0.2 

- 

- 

- 

Stellar cluster 


Comeron et 

il. 

19 

96 


ct Ori 

0.2 

- 

- 

- 

Stellar cluster 


Sherry et al. 

2( 

104) 


(4 Ori) 

(0.6) 

- 

- 

- 

Stellar cluster 


Barrado y N 

ivascues et al. 

(2004 


Note. — Listed quantities are: mass, M; the expansion velocity, v exp ; and the thermal (£ t h) and kinetic (Ekin) energy. 11 the expansion velocity in Ho and 
HI are not measured at the same position causing the discrepancy between both numbers. ®: This value has been revised downwards compared to Ochsendorf 
|& Ticlcn~s|(2015| , who determined the kinetic energy assuming an abrupt energy input appropriate for supernova remnants {Chevalicr| 1974} . Here, we argue that 
the energy input has been continuous through stellar winds (see text) and instead use E^in = l/2Mv^ xp . 


total mass within a region of r x = 160 pc would amount to M x 
= 1800 Mq. Table [3] also denotes a value for the hot gas in¬ 
side the Barnard’s Loop bubble of radius 35 pc in case it were 
filled with gas of the same character istics as observed in the 

1 .1. 2), the high column 
t prohibit 


superbubble. As noted before (Sec. 4 


densities for the Barnard’s Loop region prohibits the detection 
of diffuse X-rays towards the east side of the Barnard’s Loop 
bubble. Hot gas is detected towards the west side where the 
HI foreground column is much lower, but it is unclear if this 
emission traces hot gas from the Barnard’s Loop bubble or the 
entire Orion-Eridanus superbubble. 

Ionized gas (~10 4 K; Tablepb in the Orion-Eridanus region 
is associated with the outer Ha filaments of the superbubble, 
the Barnard’s Loop bubble shell, and the champagne flows of 
the A Ori bubble, GS206-17 +13, the IC 434 e mission nebula, 
and the Orion Nebula (Sec. |3.2| & Sec. |3.3| ). For Barnard’s 
Loop, the total ionized mass is estimated at M — 6700 Mq by 
adopting a half-sphere geometry with thic knes s 1°, outer ra¬ 
dius of 7°, and density n h = 3.4 cm -3 (Sec. 3.11. The western 
half of the complete bubble (Fig. [4| is much fainter, which 
implies that the density in the Barnar d’s Lo op bubble varies 
significantly between both parts (Sec. |4. 1.3 1 . Assuming 


= 100 km s 1 (Sec. |3.3|>, E^ m = 3.4 x 1 Cr° erg. The total 


exp 


energy that has been injected in the region can then be esti¬ 
mated at Ekin = 6.7 x 10 5() erg by assuming spherical sym¬ 
metry. This energy should be compared to the typical kinetic 
energy provided per SN (10 5l) - 10 51 erg, depending on the ef¬ 
ficiency at which it couples to the ISM; Veilleux et al. 2005 


as well as the wind mechanical energy input. By adopting 


typical wind parameters (M - 1 x 10 6 M 0 yr , v ra = 2000 
km s 1 ) for the 7 most-massive stars in the Orion OB associa¬ 
tion, and as suming the lifetime of Orion’s Cloak at 3 x 10 5 yr 
(Sec. |3.3.T]|Cowie et al.|l979| ), the total mechanical luminos¬ 
ity exerted by the Orion OB association is L mec h = 8.4 x 10 49 
erg, of which ~20% c ouples to the surrounding medium as 
kinetic energy (van Buren |T986] ). Thus, the (combined efforts 
of) stellar winds from the Orion OB 1 stars seem insufficient 
to provide the measured kinetic energy of the Barnard’s Loop 
bubble and we associate the Barnard’s Loop bubble with a 
recent supernova explosion. 

Table[3jrefiects the omnipresence of champagne flows in the 
Orion re gion, consistent with the findings of Ochsen dorf et al.| 
( |2014b) >, who showed that H II region champagne flows are 
ubiquitous throughout the Galactic plane, and illustrated the 
importance of this phase to the turbulent structure of the ISM. 
Still, the amount of material that is eroded from the molecu¬ 
lar clouds and delivered by champagne flows in the form of 
kinetic energy of the fast moving ionized gas is difficult to de¬ 
termine, and estimates over the lifetime of the star will depend 
most critically on the ionizing luminosity and main sequence 
lifetime of the ionizing star (e.g.,|Whitworth|1979)l. Here, we 
will follow |Tielens| ((2005) and adopt a ‘typical’ blister phase 
that will accelerate M G ~ 3 x 10 3 Mq of ionized gas to 30 km 
s 1 , providing a kinetic energy of in ~ 3 x 10 49 erg. We use 
this value for the champagne flow of the GS206-17+13 bubble 
that has occurred in the past (Sec. |3.3.4|i. For comparison, the 
ionized interior density of the A O ri regi on, that seems to be 
on the verge of breaking out (Sec. 333}, is quoted to be be- 
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tween 2000 M Q (van Bure n|1986j ) and 6000 M Q ( jReich|1978j ): 
this would add a similar amount of mass and kinetic energy to 
the superbubble interior as the GS206-17+13 shell did. The 
Orion nebul a has just recently entered the champagne phase 
(Sec. 3.3.4|i and only contains some ~ 20 M 0 of ionized gas 
(Wilson et al. 1997| >. The IC 434 region is a special case, as 
the champagne flow originates from a chance encounter be¬ 
tween the cr Ori star cluster and the Orion B molecular cloud 
( |Ochsendorf et al.|2014a| [Ochsendorf & Tielens|2015| >. Still, 
the star has already eroded 100 M 0 and will add another 350 
M g to the ionized gas before it plunges into the molecular 
cloud (for a current distance of 3.5 pc and a closing velocity 
of -10 km s -1 ; Ochsendorf & Tielens 2015). 

Neutra l gas (~10 z K; Tab le [3|> is present in the superbub¬ 
ble wall (Brown et al. 1995 1, and the slowly expanding shells 
of the A Ori re gion (|Lang & Masheder 11998) 1, the GS206- 
17+13 bubble (Ochsendorf & Tielens||2015| ), and th e Orion 


Nebula that is encapsulated by the neutral ‘Veil’ (e.g., van der 
|Werf et al.|20 13). The GS206-17+13 is a stellar wind-blown- 
bubble (Sec. |3.3.4| i. The Orion nebula Veil contains a total 
neutral column density of Nu = 5.8 x 10 21 cm~ 2 , has a size 
of at least 1.5 pc, as it is common to the Orion Nebula and 
M43, and shows only minor blue shifted velocities with re¬ 
spect to the background OMC (~2 km s 1 ). This implies that 
the mechanical feedback of the Trapezium through, e.g., the 
champagne flow or stellar winds has not yet coupled as ki¬ 
netic energy to the surrounding neutral gas, consistent with 
the young age of the Trapezium stars (<1 Myr; Hillenbrand] 


1997) and, perhaps, the very small dynamical time inferred 
for the Orion Nebula_since the onset of the champagne flow 
(~1.5x 10 4 yr; O’Dell etal. 2009 1 . Here, we estimate the total 


H I mass for the Orion Veil by assuming that the measured col¬ 


umn density by van der Werf et al. ( 2013] ) is homogeneously 
distributed over a shell with a radius of 2 pc, large enough 
to cover both the M43 region as well as the entire Extended 
Orion Nebula ( |Gudel et ak|2008| l. 

Molecular gas (~10 K; Table [3]) is detected in dozens of 
small cometary clouds, but largely concentrates in the Orion 
molecular clouds, which constitute the main reservoir of the 
dense gas out of which star formation occurs. As described 
by, e.g., [Bally] (|2008j), the ages and locations of the sub¬ 
groups in Orion OBI indicate that star formation has prop¬ 
agated through the ‘proto-Orion cloud’ sequentially. The ki¬ 
netic energy associated with the expansion of the A Ori molec¬ 
ular shell (6 x 10 49 erg) is consistent with expansion driven by 
either the overpressure of photo-ionized gas (Spitzer] fl978) > 
that de posits s ome 10 49 - 10 5 ° erg to the surrounding medium 
( [Ochsendorf et al!] 201 4b| ) over the lifetime of A Ori (~10 
Myr), or a S N explosion (Veilleux et al. 2005| ), or a wind¬ 
blown ori gin ( Weaver et al. |1977| , as L meC h = 2 x IQ 30 erg 
for A Ori (van Buren 1986) assuming an age of 5 Myr (Bally 
2008)1, a fraction of which will couple as kinetic energy. Some 


previous works favor a SN origin for the A Ori bubble, which 
would al so exp lain the current space motion of the A Ori star 
(see Sec. |3.3.3| l. 

Stars are listed for comparison, where we have included the 
separate subgroups from the Orion OBI association (la, lb, 
lc, and Id; |Blaauw]|1964[ |Brown et al.||1994| l. The masses 
in the la, lb , and lc subgroups have been determined using 
the results in Brown et al. ( [1994| , who constructed the initial 
mass function (IMF) and through this reported the number of 
stars initially found in the 2-120 M 0 mass range for each 
subgroup. We complement this by estimating the contribu¬ 


tion from low mas s stars below the completion limit of the 
Brown et al. ( 1994[ ) study (0.08 - 2 M 0 ) using a standard IMF 
( jKroupa 200 ip, normalized to the results from Brown et al. 
( jlPP4| ). Table [3] lists the resulting total masses, assuming an 
upper stellar mass limit set by the age of each subgroup. In 
Orion OB la, stars more massive than 13 M 0 have exploded 
as SNe, while the (younger) age of the Orion lb and lc sub¬ 
groups imply an upper mass limit of about 20 M 0 ([Bally] 
2008). Orion OB Id is the younge st su bgroup and includes 
many star-forming clusters ( |Bally|[20081 >, and here we only 
include the two largest components, the Orion Nebular Clus¬ 
ter in OMC A, and NGC 2024 in OMC B. The cr Ori group is 
noteworthy here, as it has traditionally been assigned to Orion 
OB lb based on its spatial location amongst the Belt stars, but 
its traversed space motion reveals that the cluster might have 
originated from the Orion OBlc region (Ochsendorf & Tie- 
le ns|2015"] l. We have bracketed the A Ori region in Table [3jas 
it is separated from the star clusters that are associated with 
the OMCs, and does n ot dire ctly influence the scenario that 
will be outlined in Sec. 14.3.21 

Summary: Most of the mass of the Orion-Eridanus super¬ 
bubble resides in the outer superbubble wall, while the nested 
shells that process and sweep-up the interior material (see Sec. 
4.3.2|) constitute only a small part. The molecular clouds still 
represent a large reservoir of mass, comparable to that located 
in the superbubble wall. Most of the energy is found in the 
form of kinetic energy of the expansion of the superbubble 
wall. Still, about 10 - 15% of the total energy is located in 
thermal energy of the hot gas of the superbubble interior, and 
a comparable amount lies within the kinetic energy of the ion¬ 
ized gas. The thermal energy of the hot gas created by fast 
shocks is large compared to the thermal energy of the 10 4 
K ionized gas delivered through photo-ionization from Orion 
OBI. Finally, when all subgroups in Tableware combined, 
stars comprise some ~6.9 x 10 3 M 0 or 3 - 4% of the molecu¬ 
lar material within the OMCs. 

4.3.2. Mass-loading and cleansing of the interior 

The expanding shells in the interior of the Orion-Eridanus 
superbubble are the manifestation of feedback from Orion 
OBI. When moving supersonically, the shockwaves associ¬ 
ated with the shells sweep up the interior mass of the super¬ 
bubble, acting as a mechanism to remove any material that 
may have been introduced to the medium in the past. In this 
respect, the mass (6.7 x 10 3 M 0 ) and radius (35 pc) of the 
shell surrounding the Barnard’s Loop bubble can be combined 
to show that the average density of the swept-up medium con¬ 
tained ~1 cm -3 . We note that the density of the swept up 
material only slightly exceeds what would be expected for 
the WNM of the Galaxy, which would be consistent with the 
Barnard’s Loop bubble being a superbubble wall expanding 
into the ambient Galactic medium. However, we render this 
scenario unlikely given the short dynamical timescale of the 


Barnard’s Loop bubble (-0.3 Myr; Sec. |3.3.2|i compared to 


the Orion star forming region (-15 Myr; |Ba 
pre-shock density of the HV gas (~ 3 x 10~ J cm 


ly 2008), the 
— -jt 


Vfelty 


et al. 20021 and its nested appearance within the entire Orion- 


Eridanus superbubble (Figs. [I]&|4|. The increased density 
is about two orders of magnitude above that what is expected 
for a superbubble interior (~ 10 2 cm~ 3 ; Weaver et al. 1977), 


which provides convincing evidence that the superbubble in¬ 
terior must have been mass-loaded since the last supernova 
passed through and cleansed the interior. 

Photo-ionization feedback from Orion OB 1 ultimately ere- 
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ates champagne flows that channel 10 4 K gas from the Hll 
regions into the superbubble cavity, where it will mix with 
the hot gas (>10 6 ; Chu 2008 ). The injected mass will signif¬ 
icantly lower the temperature of the superbubble interior and 
increase its density. The total mass delivered to the interior by 
a typical champagne flow over the lifetime of a massive star 
(~3 x 10 3 M 0 ) only slightly exceeds the mass that makes up 
the hot gas of the entire superbubble (1.8 x 10 3 M G ). How¬ 
ever, this extra mass will be injected in a small portion of 
the interior (see below). In addition, champagne flows will 
dump kinetic energy in the form of turbulent motion into the 
hot cavity. In principle, the injected turbulent energy is more 
than the thermal energy of the hot gas provided by the SN that 
lead to the formation of the Barnard’s Loop bubble (Table[3]t. 
However, dissipation of this turbulence will lead to tempera¬ 
tures of only ~ 10 4 K (for a flow speed of ~ 30 km s _1 ; Table 
[3j Tielens ]|2005ll. The cooling function peaks at these tem- 
peratures (|Dalgamo & McCray 1972) such that this excess 


energy wilt be mainly radiated away. Nevertheless, adding 
and mixing in the champagne gas will mass-load and cool the 
hot superbubble interior gas. In addition, thermal evaporation 
from the dozens of clouds that are submersed in the hot gas 
of the superbubble, most notably the Orion Molecular clouds, 
will contribute to the mass-loading of the hot gas. It is diffi¬ 
cult to constrain the relative contributions to the mass-loading 
of the hot gas through thermal and photo-evaporation, as the 
therm al evaporation rate is a s teep function on temperature (oc 
T 5 ' 2 - Cowie & McKee 1977), while photo-evaporation rates 
depend on the ionizing luminosity a nd main sequence lifetime 
of the ionizing star (e.g., Whitworth|1979| 

er that thi 


1979 ). Nonetheless, we can infer 


Bodenheimer et al. 
e combined effects 
of the mass-loading mechanisms have introduced some 6.7 x 
10 3 M 0 of material into the superbubble cavity over the past 
3 x 10 5 yr given by the mass and age of the Her shell associ¬ 
ated with the Barnard’s Loop bubble that acts to cleanse the 
interior of the superbubble close to the OMCs. These num¬ 
bers imply an average evaporation rate (or cloud destruction 
rate) of ~1 x 10 2 M 0 yr -1 , which is an order of magnitude 
above the current photo-erosion rate in the Orion Nebula, if 
we use a total ionized gas mass of 20 M 0 (Wilson et al. 1997) 
and a dynamical lifetime of 1.5 x 10 4 yr (O’Dell et al. 2009 1 . 
Part of this discrepancy may be because the Orion nebula still 
resides in an early phase of evolution and the region has not 
fully broken out of its natal cloud. Nonetheless, mass-loading 
of the superbubble interior increased the density by a factor of 
~100 compared to that expected in case the region within the 
Barnard’s Loop bubble were filled with hot gas of the same 
characteristics as observed in the superbubble (TableRl). Over 
time, this will have a rigorous effect on the thermal behavior 
of the interior gas. 

We note that the pre-shock density derived for the HV 
gas associated with the Ba rnard ’s Loop bubble in the line of 
sight towards ( Ori (Sec. 3.3 1 is 3 x 1CL 3 ctrT 3 , while we 
have derived an average density of 1 cm -3 for the medium in 
which the Barnard’s Loop bubble has expanded in its 3 x 10 5 
yr lifetime. Surely, the mass-loading mechanisms described 
here concentrate towards the centre of star formation, i.e., the 
OMCs. There are believed to have been 10 - 20 SN explo¬ 
sions within Orion in the last 12 Myr ( ]Bally|j2008j ). With 
an expansion velocity of -30 km s 1 for the gas introduced 
through photoablation and champagne flows (Table [3}, a su¬ 
pernova rate of 1 - 1.5 Myr -1 will not allow the injected mass 
to disperse and fill the entire volume of the superbubble before 


being swept-up by a subsequent SNR. Thus, the increased 
densities due to mass-loading from the OMCs will be mainly 
concentrated within 30 - 45 pc surrounding the OB associa¬ 
tion (consistent with the radius of the Barnard’s Loop bubble), 
whereas the surrounding medium will have a much lower den¬ 
sity. The low temperature in the eastern part of th e supe rbub- 
ble between Barnard’s Loop and region 1 (Sec. |4.1.3|) may 
also be caused by the effects of mass-loading. This would im¬ 
ply that other ‘routes’ of mass-loading affect the density and 
temperature structure of this region, such as enclosed clouds 
toward the Galactic plane, the straight evaporation of the su¬ 
perbubble wall, or SNRs that dissipate before reaching the 
superbubble wall (see below). However, as we have not been 
able to constrain the density structure of this region, the im¬ 
portance of mass-loading remains unclear. 

In contrast to the shell of the Barnard’s Loop b ubble that 
is kept photo-ionized by Orion OBI (Sec. |3 .2.2)), the neu¬ 
tral shells surrounding A Ori and GS206-17+13 have had the 
chance to cool radiatively and collapse to form dense media. 
While they may add to the mass-loading of the superbub¬ 
ble interior through thermal evaporation, their slow expansion 
renders them unimportant compared to the total energetics of 
the superbubble (Table[3]t. Still, these dense shells may be the 
site for a new generation of stars a s they are prone to gra vita¬ 
tional instabilities (Zavagno et al.|2007| Pomares et al.|20d9| 


|Martins et al.|2010'. 

The supernova remnant associated with the Barnard’s Loop 
bubble has exploded ~3 x 10 5 yr ago and currently expands 
at 100 km s , implying that it is well in the momentum- 


conserving snow plow phase (e.g., 
shock gas contains densities (n — 3. z 

Tielens 

2005). 

The post- 
) and tem- 

cnT 3 ; 

3ec.[37l 

peratures (T = 6000 K; Heiles et al. 

2000 O’DeTTet al.|2011 


Madsen et al.|2006) that would allow the gas to cool rapidly 


(Tcooi ~ kTl nA(T) ~ 4 x_10 4 yr, where A (T) » 6 x 10 
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erg cm 3 s 1 ; 


Dalgamo & McCray 1972). However, Barnard 


Loop is kept photo-ionized by Orion OB 1 (Sec. 3.2.2 1 , pro¬ 


hibiting the cooling and collapse of the gas to form a dense 
neutral shell. In this case, the shockwave will continue to 
gather mass until it becomes subsonic, from which point it 
will dissipate and pressurize the interior, thus effectively dis¬ 
placing material from the origin of the blast wave towards the 
outskirts of the superbubble interior. 

To get a feeling how far the displacement of material will 
reach, we can compare radii and expansion velocities of the 
Barnard’s Loop bubble (r = 35 pc, v exp = 100 km s 1 ) and the 
superbubble (r = 160 pc, v exp = 40 km s" 1 ) to estimate that 
the Barnard’s Loop bubble will catch up with the outer wall 
within ~2 Myr if both structures decelerate at the same rate 
(in the radiative expansion phase, the velocity of the SNR de¬ 
celerates with time as oc t 3/7 ; Tielens 2005). At this moment, 


the isothermal sound speed in the outer regions of the super¬ 
bubble towards the east is ~10 km s -1 (for T ~ 10 4 K; Sec. 
|4.1.3| >. If one assumes the lifetime of Barnard’s Loop at 3 x 
10 3 yr, and adopt the time-velocity relation of the SNR above, 
the adiabatic expansion phase (v exp > 250 km s 1 ) ended just 
below 10 s yr and it will take another ~7 Myr to reach 10 km 
s if the density ahead of the SNR is homeogeneous. This 
timescale far exceeds that derived for the Barnard’s Loop bub¬ 
ble to catch up with the superbubble wall. Thus, Barnard’s 
Loop will be able to cross the superbubble radius and reach 
the outer wall, where the shockwave will eventually thermal- 
ize and the swept-up material will transfer its mass and mo¬ 
mentum by condensing onto the superbubble wall (Mac Low 
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& McCray|jl988)). The outer superbubble wall may therefore 
be accelerated as various interior shells such as Barnard’s loop 
catch up to it. In between such pulses, it could be coasting and 
slowing down as it sweeps up surrounding gas in a momen¬ 
tum conserving interactio n, or i t may be dissolving in case it 
moves subsonically (Sec. |4.1.3|). 

We note that densities in the outer parts of the superbubble 
are smaller compared to the mass-loaded region close to the 
OB assocation (~ 1 cm' 3 ; see above). In this case, the result¬ 
ing pressure drop will accelerate the interior SNR once again 
when it reaches the edge of the mass-loading region, which 
already seems to be the case for the Barnard’s Loop bubble, 
given the pre-shock density of 3 x 10 3 cm' 3 in the line of 
sight towards £ Ori ( [Welty et al.|2002 i. Indeed, towards the 
west side of the superbubble, the temperature of the interior 
is much higher and the SNR will accelerate, but it may go 
subsonic and dissolve earlier compared to the east side of the 
bubble: this may be at the root o f the fa intness of the Her emit¬ 
ting shell in this direction (Sec. |4.1.3|). Towards the east, the 
density gradient could be smaller because of the lower tem¬ 
perature that may be accompanied by an increase in density, 
and the amount of acceleration will decrease accordingly. 

If 10 - 20 SN explosions occurred in the Orion region in the 
past 12 Myr ( Bally|20 08 ), and each SNR adds ~7 x 10 3 Mg 
of material to the superbubble wall (Table|3j, this will amount 
to 0.7 - 1.4 X 10 5 M g and account for a significant fraction 
(20 - 40 %) of its total mass. Thus, the mass of the super¬ 
bubble wall without the addition from SNRs would account 
for 1.9 - 2.6 x 10 5 M 0 . With a radius of 160 pc, we can esti¬ 
mate the average ISM density into which the superbubble has 
expanded at n ~ 0.15 - 0.2 cm' 3 , consistent with the WNM 
within one or two scale heights of the Galactic plane. Finally, 
as the supernova rate is of order for the Barnard’s Loop bub¬ 
ble to catch up with the outer wall, it may not be necessary to 
expect multiple layers of expanding shells between the Orion 
OB association and the superbubble wall. 


4.3.3. Efficiency of star formation in the Orion molecular clouds 


Given a total mass of ~2 x 10 5 M- for the OMCs combined, 
a supernova rate of 1 - 1.5 Myr -1 (Bally 2008) and swept-up 
mass of ~1 x 10 3 Mq per SNR (the mass of the Her emitting 
part of the Barnard’s Loop bubble), we estimate that at this 
rate the reservoir of cloud material will be depleted within 20 
- 30 Myr through thermal evaporation and erosion from cham¬ 
pagne flows alone, similar to calculated cloud lifetimes from 


Williams & McKee (1997). Thus, material is effectively re¬ 
moved from the molecular clouds and incorporated in expand¬ 
ing shells that might form stars of their own. Furthermore, the 
molecular clouds will be disrupted by each supernova blast 
wave that mark the end of a cycle of star formation. All of 
these effects should be considered if one attempts to simulate 
the lifecycle of molecular clouds and the efficiency of star for¬ 
mation, which is well known to be limited by various forms 
of stellar feedback (for a recent review, see IKrumholz et al.l 
[20l4|. 


The gradual erosion of the OMCs puts a firm upper limit 
on the lifetime of the molecular clouds. Over the past 10 - 15 
Myr, some ~10 5 Mq has already been lost from the OMCs to 
the superbubble wall, which implies that the original molec¬ 
ular reservoir (the ‘proto-Orion’ cloud) contained some ~3 X 
10 5 M 0 (Table a At present, the major subgroups_of the 
Orion OB association (B laau wj 1 964| |Brown et al.|1994| ) con¬ 
stitute some 6.9 x 10 3 M Q of stellar mass. The ages of the 


subgroups impl y tha t many stars have already exploded as su¬ 
pernova (Sec. |4.3| |Bally||2008j), and we estimate from the 
adopted IMF that about 20% of the initial mass of the la, lb, 
and lc subgroups has been lost through SN explosions, in¬ 
creasing the total mass that has been incorporated into stars 
over the last 10 - 15 Myr towards ~ 7.4 x 10 3 M 0 . With a 
remaining cloud lifetime of the OMCs of 20 - 30 Myr, at the 
current rate, some ~2.2 x 10 4 M 0 or 7.5% of the molecu¬ 
lar mass can be converted into stars over the lifetime of the 
clouds con sistent with predictions by, e.g., [Williams & Mc-| 
Kee|(|1997j). However, it is unclear if star formation will pro¬ 
ceed similarly as stellar feedback may trigger star formation 
and disrupt the clouds, while the molecular reservoir steadily 
diminishes (see, e.g., Evans et al.|2009. Hopkins et al.|201 1| 


Krumholz et ak 2014 for more thorough discussions regarding 

the mechanisms that involve star formation rates and molecu¬ 
lar cloud lifetimes). 


4.3.4. Summary: the rejuvenation of the Orion-Eridanus 
superbubble 

We have argued that appreciable amounts of mass will be 
removed from the molecular clouds and be mass-loaded into 
the interior of the superbubble through destructive champagne 
flows and thermal evaporation of clouds embedded in the hot 
gas of the superbubble interior. Explosive feedback from stel¬ 
lar winds and, in particular, supernovae take care of the subse¬ 
quent transportation of the mass-loaded material towards the 
outskirts of the superbubble through the formation of expand¬ 
ing nested shells. For the Barnard’s Loop bubble towards the 
east, photo-ionization by Orion OBI does not allow it to radia- 
tively cool and collapse, and it will plaster the mass and add 
momentum to the superbubble wall that will continue to drive 
the expansion of the superbubble. In case it would go sub¬ 
sonic before reaching the outer wall, the shell would dissolve, 
thermalize, and pressurize the interior. This cycle repeats as 
long as there is ongoing star formation and OB stars can main¬ 
tain (the photo-ionization of) the cleansing shells that deposit 
the mass towards the outer regions. 


4.4. Dust processing in superbubbles 

An interesting consequence of the constant mass-loading 
of the superbubble is that each supernova explosion within 
the superbubble cavity will interact with a significant amount 
of ‘fresh’ interstellar gas that carries along dust in its wake 
through entrainment. The efficiency of dust entrainment in 
photo-evaporation flows was investigated by Ochsendorf & 
Tielens (2015) for the IC 434 champagne flow, revealing that 
the dust-to-gas ratio within the flow is similar to that seen in 
the diffuse ISM (~0.01). Thus, each single SNR from the 
Orion OBI association will encounter a new reservoir of gas 
and dust within the superbubble cavity, introduced by cham¬ 
pagne flows from the HII regions and from photo-evaporation 
flows of the thermally evaporating clouds. This finding con¬ 
trasts with that of|McKee|(|1989|), who previously estimated 
the timescale for SNR sh ocks to destroy a typical dust grain 
within the ISM. McKee ( 1989| > assumed that only the first su¬ 


pernova within an OB association would destroy an amount of 
dust equivalent to that located in an ISM volume containing 
~1300 M 0 of gas. Following supernovae from the OB associ¬ 
ation would hammer surrounding gas already cleansed of its 
dust component. Following [McKee] (|l989j), the effective SN 
rate is equal to 
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‘SN 


= (?l/l + [/lI,0B<7lI,0B + /lI,field?II,field] /ll)^SN yT. (6) 


Here, <71 is the correlation factor of the SNe of type la 
(~0.38; Heiles|1987[ i, relating the location of the SN to their 
position above the disk of the Galaxy (i.e., it measures the ef¬ 
fective interaction of the SNe with the amount of gas in their 
environment^ and fi is the fraction of supernovae of type la 
(~0.5; Narayan |l987) >. Similarly, fi\ (= 1 - fi) measures the 
fraction of supernova of type II, part of w hich occur in OB as- 
sociations, described by /h,ob (~ 0.75; see Zinnecker & Yorke 
|2007[ and references therein), and in the field, denoted with 
/n.fieid (= 1 - /ii.obX originating from runaways and isolated 
stars with correlation factor < 7 n,fi e id (— 0 . 6 ; using values for the 
scale height of the WNM by Heiles|1987) . Finall y, *~cm is the 


Diehl et al. 


intrinsic SN rate in the Galaxy (~2 x 10 " yr 
|2006| . 

Under the assumption that superbubbles do not replenish 
their dust content after the first SNR, |McKee| ( |1989| arrived at 
a correlation factor of c/h ob = 0.1 and, together with a rather 
low value for the fraction of type II SNe in OB associations, 
/ii.ob = 0.5, derived an effective SN rate of tsn = 125 yr. 
In this work, we have shown that dust entrainment accom¬ 
panying mass-loading of the superbubble interior replenishes 
the dust content inside the superbubble after each SNR from 
the OB association, i.e., < 7 n,OB = 1. In addition, we use the 
more recently derived fraction of OB stars that are located 
in OB associations from Zinnecker & Yorke (2007), /h ob = 
0.75. These changes lead to shorter effective supernova rate of 
tsn = 90 yr, implying that the SNRs within superbubbles de¬ 
stroy dust more efficiently and, subsequently, the dust lifetime 
against SNR shocks must be revised downwards by 30%. 

The processing of dust inside the superbubble can also be 
probed directly from depletion studies. |Welty et al.| ( |2002) 
noted that in the HV and IV gas towards £ Ori, the gas-phase 
abundances of Al, Si and Fe are slightly elevated compared 
to that observed in warm, diffuse clouds. However, the car¬ 
bon abundance is significantly raised in the HV gas to near¬ 
solar abundances. For comparison, depletion studies of the 
Cold Neutral Medium (CNM) and the WNM show a pat¬ 
tern where a thin coating on interstellar dust gets sputtered 
and re-accreted, as the grains cycle from clouds to the diffuse 
phase of the ISM ( |Tielens|1998]|20T3) . Whilst in the WNM, 
the reservoir of gas and dust typically encounters a 100 km 
s 1 passing shock, destroying some 10 - 30% of the silicate 
volume (Jones et al .|1996| l. However, these shock velocities 
would only destroy 15% of the carbonaceous dust and lead 
to minor variations in the carbon depletion, as ~ 50% of the 
carbon is in the gas phase (|Cardelli et al.||1996[ |Sofia et al. 
2004). Therefore, instead of passing shocks, we attribute the 


enhanced carbon abundance to sputtering of the grains inside 
the hot gas of the Orion-Eridanus superbubble. The sputtering 
timescale T sput (Tielens et al | 1994j i in the X-ray emitting gas 
(T ~ 10 6 K, n ~ 0.01 cm *) implies that large 0.1 pm silicate 
grains are hardly affected (T sput > 20 Myr), yet, small ~ 10 A 
(carbonaceous) grains are quickly eroded (r sput < 1 Myr), re¬ 
leasing the carbon into the gas-phase. We argue that this may 
be a general characteristic for dust that is processed inside su¬ 
perbubbles. Indeed, dust entering the heliosphere also shows 
a significant increase of gas-phase carbon, implying that the 
grains are sputtered within the hot gas of the Local Bubble 
( Frisch & Slavin|2013| l. 


5. SUMMARY 

The general picture that emerges from the multitude of data 
is that the Orion-Eridanus superbubble is larger and more 
complex than previously thought, consisting of a series of 
nested shells which we depict in Fig. [ 6 ] With the data cur¬ 
rently at hand, we have discussed the following structures, 
ordered in size from bottom up: 


- The youngest shell surrounds the Orion Nebula cluster 
located in Orion A (age < 1 Myr; |Bally|2008 1 . 

- Next, there is the 2 to 4 degree HI and dust shell GS206- 
17+13, approximately centered on the cr Ori cluster, but 
likely blown by stellar winds from the Belt stars (age ~ 
5 Myr; |Ochsendorf & Tielens|2015] l. 


Barnard’s Loop is part of a complete bubble structure 
not related to the outer edge of the Orion-Eridanus su¬ 
perbubble. It is instead associated with a supernova 
remnant (age ~3x 10 s yr), and is connected with high- 
velocity g as (‘ Orion’s Cloak’) detected in absorption 
studies (Cowie et a k||T979l [Weity et al.|2002) . All of 
the previous shells are projected in the interior of this 
young 7 degree radius bubble, which expands at high 
velocity (v exp ~ 100 km s _1 ). 


The Barnard’s Loop bubble and the A Ori bubble (~ 5 
Myr; |Bally|2008j) expand in the interior of the ~ 45 x 45 
degree shell that is related to intermediate-velocity gas 
from the aforementioned absorption studies and created 
by the collective effects of Orion’s stars: the Orion- 
Eridanus superbubble (age ~ 5 - 10 Myr). 


As we can recognize several distinct subgroups and well 
defined clusters, star formation in the Orion-Eridanus region 
must have been highly episodic over the last 10-15 Myr. It 
may therefore not be surprising that the Orion-Eridanus super¬ 
bubble consists of a set of nested shells. Each ‘burst’ of star 
formation may have produced a subgroup consisting of sev¬ 
eral clusters that disrupted parts of the pre-existing molecular 
cloud. The formation of each subgroup and the resulting UV 
radiation can ionize large amounts from the dense reservoir of 
molecular gas in the region that, together with thermal evapo¬ 
ration of the molecular clouds, constantly mass-load and cool 
the hot gas inside the superbubble cavity. Stellar winds and 
supernovae accelerate, sweep-up, and compress these ‘poi¬ 
soned’ plasmas in an episodic fashion to form nested shells 
within the Orion-Eridanus superbubble, such as the Barnard’s 
Loop bubble, the A Ori bubble, the GS206-17+13 shell, and 
the Orion nebula (see Fig. [ 6 ]). The shells may cool, collapse, 
fragment, and be incorporated in a next generation of stars. 
However, for the Barnard’s Loop bubble, this shell-cooling 
is inhibited by photo-ionization from the OB association. We 
have noted that in its turn, the Orion-Eridanus superbubble has 
previously been associated with the ancient supershell known 
as Lindblad’s ring, which may have cooled and collapsed to 
form the Gould’s Belt of stars that includes the Orion OB as¬ 
sociation, extending the hierarchy of nested shells and bubbles 
towards greater sizes and into history. The continuous replen¬ 
ishment of dust in the superbubble cavity through entrainment 
in ionized flows has lead us to conclude that dust process¬ 
ing from interior supernova remnants is more efficient than 
previously thought. To conclude, the cycle of mass-loading. 
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cleansing, and star formation ceases when feedback has dis¬ 
rupted the molecular reservoir, from which the superbubble 
will disappear and merge with the ISM. 
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